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TRANSITION ASYMPTOTICS FOR THE PAINLEVE II TRANSCENDENT 

THOMAS BOTHNER 


Abstract. We consider real-valued solutions u = u{x\s),x S M of the second Painleve equation u^x = 
XU + 2u^ which are parametrized in terms of the monodromy data s = (si, S2, S3) C of the associated 
Flaschka-Newell system of rational differential equations. Our analysis describes the transition, as a: —> — 00 , 
between the oscillatory power-like decay asymptotics for |si| < 1 (Ablowitz-Segur) to the power-like growth 
behavior for |si| = 1 (Hastings-McLeod) and from the latter to the singular oscillatory power-like growth for 
I SI I > 1 (Kapaev). It is shown that the transition asymptotics are of Boutroux type, i.e. they are expressed 
in terms of Jacobi elliptic functions. As applications of our results we obtain asymptotics for the Airy kernel 
determinant det(/ — 7 RAi)li, 2 (j. in a double scaling limit x —> — 00,7 f 1 as well as asymptotics for the 
spectrum of i^Ai. 


1. Introduction and statement of results 


One of the most impressive occurrences of Painleve transcendents in non-linear mathematical physics stems 
from their applicability in random matrix theory. In this field Painleve functions describe, for instance (cf. 
[32, 3, 18]), eigenvalue distribution functions for classical finite n ensembles, they appear in the description of 
universal distribution functions in the large n limit and are also used in the computation of gap probabilities 
in the large n limit. One concrete example for the last two cases is given by the celebrated Tracy-Widom 
distribution [31]: the distribution function for the largest eigenvalue A^ax of a n x n random matrix drawn 
from the Gaussian Unitary Ensemble (GUE) in the large n limit equals 


Ftw(t) = lim Prob ( Amax < + 

lastin 

0 = 


= exp 


iy - x)ul^{y)dy 


X G 


n-i-oo ^/2n^ 

where Uhm = Uhm{x),x G M is the Hastings-McLeod solution of the second Painleve equation, 

d , ^ x~i 


dx’ """" (1 +'^(1))’ 


-boo. 


( 1 . 1 ) 


( 1 . 2 ) 


This special solution was first analyzed in the works of Hastings and McLeod [21] who showed that the 
boundary value problem (1.2) has a unique, monotonically decreasing, smooth solution. In addition. 


,{x) 


x^- 00 . 


Although Ftw(^c) in (1-1) is written as a distribution function, it is equivalent to a gap probability for 
the rescaled eigenvalues ftj = •\/2r®(Aj — and as such, expressible as Fredholm determinant. More 

precisely, 

FTw(a::) = J.^Prob(^|l{j : G [x,oo)} = oj = det (I - Aai) , (1.3) 

where Aai is the trace-class operator in L^((a;, 00 ); dA) with kernel 


(XjOo) 


Ai(A)Ai'(bf) - Ai'(A)Ai(/r) 

AAi(A, /r) = ' G A, G (x, 00 ) 

A — fl 
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and Ai(A) the classical Airy function. Besides the gap probability (1.3) other spectral properties of large 
random matrices from the GUE can also be computed through Painleve transcendents, for instance for 
k G Z>o we have (cf. [28]) 


lim Prob 

n—>-<50 


(i{j : G [x,oo)} = fc) = ^ r- 


^7 


k r 


det (/ - yATAi) 


(ai,oo) 


7=1 


and 


det (/ - yATAi) 


= exp 


{y - x)ulsiy;-f)dy 


y > 0. 


(rc,oo) 

Here, Uas = Uas{x]j),x G K is the Ablowitz-Segur solution of the second Painleve equation. 


m" = XU AS + 2ui 




UAsix;^) = 7y 


2^/n 


9 ^ 

e-§- = 


(l + o(l)), X +00. 


(1.4) 


(1.5) 


This one-parameter family of solutions has different analytical and asymptotical properties depending on 
the values of y: namely, for the values of y G (0,1) fixed, the boundary value problem (1.5) has a unique, 
bounded, smooth solution with oscillatory behavior as x —>■ —oo, compare (1.18) below. On the other hand, 
if y > 1, smoothness is destroyed at hnite x and the solution blows up, see (1.20). 

Remark 1.1. Historically, solutions of the boundary value problem (1.5) for fixed y G (0,1) were first 
analyzed, both analytically and asymptotically, in the late 1970’s, cf. [1, 2]. In particular for the bounded so¬ 
lutions, Ablowitz and Segur solved the important connection problem, i.e. determine the complete asymptotic 
description of Uas(x]^) as x ^ —oo provided the same description is given as x ^ -|-oo (or vice versa). The 
same problem was subsequently also solved for the unbounded solution Mhm)^;) of (1.2), see [21]. However, 
the singular asymptotic structure for y > 1 as x ^ —oo remained unknown until the work of Kapaev [26] in 
1992. 

Our modest goal in this paper is to place Ablowitz-Segur (for y G (0,1) or y G (l,oo)) and Hastings- 
McLeod solutions on equal asymptotic footing as x —> — oo. This problem is of interest to the asymptotics 
of the Fredholm determinant 

det(/-yA:Ai) 

(ic,oo) 

as X —>■ —oo and thus to the asymptotics of the Tracy-Widom distribution itself. 

For the gap probability let us use the refined asymptotic behavior (compare Remark 1.8 below). 




( 1 . 6 ) 


and the identity 


- u„ 


POO 

/ Kmiy)dy = - xul^ix) 

J X 

from which we learn in (1.1) that 

-®co(l + 0(x-3)). 


(a;), x G 


FTw(a:) = exp 


12 


(1.7) 


with some constant Cq G K. The constant factor was first derived by Deift, Its and Krasovsky in [11] using 
scaling limits for a finite Laguerre ensemble combined with universality results for the Airy kernel, 


Co = exp 


^In2 + C'(-1) 


( 1 . 8 ) 


with C = f{.s) the Riemann-zeta function. Shortly afterwards Baik, Buckingham and DiFranco [4] gave 
another derivation of (1.7), (1.8) based on exact integral representations of FTw(a;) which, opposed to 
(1.1), involve an integration from —oo to x. In either approach the regularization of the integral in (1.1) 
is important and thus the asymptotics (1.6) relevant as they are needed in the evaluation of a large gap 
probability in random matrix theory.^ 


^Another way of deriving (1.7) without the Painleve II connection was presented in [6]. 
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The very same asymptotic question for the behavior of det(/ — yi^Ai) with fixed 7 7 ^ 1 is still open and 
our focus does not lie on this interesting problem. However, considering the different qualitative behaviors 
of Uj^s{x;j) as X ^ —00 (see (1.18), (1.20) below for further detail), 


Mas (a;; 7) 


i-x)i 


COS 


+ /31n(8(-a;)2 )+0 (^(-x) 10 ^ 


3(-a^) = 


= ^ln(l- 7 ), </>= ^-argr(i/3); 


ye (0,1) 


Mas (a:; 7 ) 


sin (|(—x )2 + /31n(8(—x) 2 ) + (/?)+ C>((—x) 2 ^ 


+ 0((-x)-i). 


7 G (1,00) 


^=^ln( 7 -l), tp = ^ - argT Q+i^^ ; 

it becomes evident that det(/ — 7 ^fAi)|L 2 ( 2 , experiences a phase transition near 7 = 1 as x —>■ — 00 : for 
7 G (0,1) the determinant will be strictly positive, decaying exponentially fast but with a slower decay 
rate than (1.7). Opposed to that, for 7 G (l,oo) the determinant will display oscillations with decreasing 
amplitudes and its zeros are accumulating at x = — 00 . This qualitative change leads to a challenging 
problem: Determine the asymptotic behavior of 


exp 


(y - x)ul^{y;'r)dy 


det (J-yiGAi) 


(rr,oo) 


as X —>■ — 00 and simultaneously 7 —)■ 1. (1.9) 


As a direct application of Theorem 1.12 below on the transition asymptotics of Uas(x; 7 ) we will prove a short 
result which shows that (1.7) to leading order also describes one particular case of the transition asymptotics 
for det (I - jKm). 


Corollary 1.2. As x 


—00 and 7 t 1 ) with cq as in ( 1 . 8 ), 


det (/ - yATAi) 


L‘^{x,oo) 


= exp 


12 


|x| sco(l + o(l)). 


uniformly for 


ln(l - 7 ) 

/ \ 3 


> 



For the values of x and 7 such that 0 < x: < our results in Theorem 1.10 and 1.13 do in general not 
allow us to state an analogous expansion; the integration in (1.9) simply becomes too complicated. Still we 
can make a qualitative prediction 

The transition asymptotics of det(/ — 7 ArAi)|L 2 (a;,oo) as x —>■ — 00,7 f 1 and 0 < x: < is 
modeled by quasi-periodic functions, to leading order. 


Remark 1.3. We would like to point out that an asymptotic analysis similar to (1.9) was recently carried 
out for the case of another prominent Fredholm determinant in random matrix theory, 


P(s, 7 ) = det(/ - yATsin) 




with Ksin{X,y) 


sin s(A — y) 
7 r(A- /r) 


s ,7 > 0 . 


It is well known [28, 32] that P{s, 1) equals a gap probability: the probability of finding no eigenvalues in the 
interval (— 7 , 7 ) for a random matrix chosen from the GUE, in the bulk scaling limit with mean spacing one. 
Similar to the Airy kernel determinant, P{s,^) admits a representation in terms of Painleve functions, this 
time involving a special solution of the fifth Painleve equation [25]. Building on previous work of Dyson [15], 
the analysis of P{s,j) as s ^ + 00,7 f 1, was essentially completed in [7], however without using Painleve 
asymptotic analysis. This lead to results completely analogous to Corollary 1.2 and the qualitative prediction, 
compare Theorems 1.4 and 1.12 as well as Corollary 1.13 in [7]. 


We now move ahead and present the setup for our results on the transition asymptotics as x —>■ —00 
between Uas{x;^) and MHM(a:). Opposed to introducing Painleve II functions as solutions to boundary value 
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problems as in (1.2) and (1.5), we choose to follow a Riemann-Hilbert point of view. In this approach [17, 24] 
solutions u = u(x) to 

Uxx=XU + 2u^ ( 1 - 10 ) 

are parametrized through the monodromy data of an associated linear system of ordinary differential equa¬ 
tions in the complex plane. The details are as follows, compare [19]: Given six complex numbers {s/c}fc=i 
which satisfy the relations 


Si - S 2 + S 3 + S1S2S3 = 0 , 
we introduce the triangular matrices 


^k-\-3 — 


Si = 53 , S2 = 52 , 


( 1 . 11 ) 


S'*: = ( ^ j for k = 1 mod 2, 


Sk = 


1 Sk 
0 1 


for k = 0 mod 2, 


and the sectors 


= |a S C : —(2fc — 3) < argA < —(2fc — 1)|, A: = l,...,6 
with oriented boundary rays 

rfc = |AeC: argA = ^-h ^(fc - 1)|, k = l,...,6 
as shown in Figure 1 below. Consider now the following Riemann-Hilbert problem (RHP). 


Riemann-Hilbert Problem 1.4. Determine the piecewise analytic 2x2 matrix-valued function Y(X) = 
Y (A; x,s = (si, S2, S3)) such that 

• Y(A) is analytic for A G C\M?rfc and for every k, the function Yfc(A) = y(A)|„ has a continuous 
extension on the closure Vlk- 


(; ■=■;”) 



1 0 

s^e^o 1 


1 See 
0 1 


1 0 
see^^ 1 


Figure 1. Jumps in the Painleve H Riemann-Hilbert problem 


• The boundary values F+(A) (resp. Y-{\)) from the left (resp. right) side of the oriented contour Tk 
are related via the jump condition 

y+(A) = A G Pfe (1.12) 

where 

6»(A,a;) = i Qa^- hxA^ , ^3 = • 

• Ts A tends to infinity, the function Y (A) is normalized as follows, 

Y{X) = I + 0{X-^), X-^ 00 . 


(1.13) 
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Remark 1.5. The first two of the relations in (1-11) ensure consistency of the jump conditions (1-12) with 
the continuity ofYk{X) on fife. Indeed, conditions (1-12) can be rewritten as 

n+i(A) A G fc = l,...,5 

ri(A) = A G Tg 

and thus 

Yi{X) = S 1 S 2 ■ ■ ■ ■ ■ , A G Ti. (1.14) 

But every solution Y[X) of the RHP l.f is invertible since detF(A) is an entire function normalized to unity 
at infinity. Henee (1-14) yields upon cancellation o/Yi(A) the cyclic constraint 

S,S 2 ■...■Se = I (1.15) 

which is equivalent to the first two conditions in (1.11). We also mention that by triangularity of Sk, we 
have ^ I as A —>■ (X) along the rays Ffc, i.e. (1-12) is consistent with (1.13). 


The connection between the RHP 1.4 and (1-10) is as follows. It is known [5] that for any set of parameters 
s = (si, 52 , 33 ) satisfying (1.11) the RHP for H(A) is meromorphically (with respect to x) solvable. Moreover 
its solution determines the Painleve H transcendent via 


u{x) 


= M(a:|s) = 2 lim 

A-s-oo 


A(y(A;x,5))^ 


(1.16) 


and in addition we have u{x) = u{x). Conversely every real for real x solution of (1.10) has a unique 
Riemann-Hilbert representation (1.16) for suitable s and we therefore adopt the notation u{x) = u(a;|s), 
indicating the parametrization of solutions to PH equation (1.10) in terms of the data s. 


Remark 1.6. The previously considered cases (1.5) and (1.2) are special cases of (1.16), 
Uas(x; 7 ) = u(x|(-i^/ 7 , 0 , 1 ^ 7 )), Uhm(x) = u(a:|(-i, 0 ,i)). 


The data s is equal to the monodromy data of the associated linear 2x2 matrix ODE 


d^' 

dA 


—diA^CTs + 4iA 


0 u 

—u 0 


—ia; — 2iu^ 

2u,x 


— 2 ux 
c + 2iu^ 


4- 


(1.17) 


where u satisfies (1.10). In the language of classical monodromy theory of differential equations (cf. [22]) the 
entire functions 'I'fe(A) = Ffc(A)e“®*-'''’®^'^ 3 ^ fc = 1,..., 6 and 4 ' 7 (A) = 'I'i(A) form the seven canonical solutions 
of the ODE (1.17) and the triangular matrices Sk in (1-12) are the Stokes matrices characterizing the Stokes 
phenomenon of the irregular singular point A = 00 . The matrices Sk are a:-independent and thus the second 
Painleve transcendent u = M(a;|s) describes the isomonodromy deformations of system (1.17) with respect 
to the (isomonodromic) deformation parameter x. 


As we are interested in the asymptotic behavior of real solutions m(x|s) as x —>■ — 00 , it will be useful to 
recall the following three different cases which have been singled out over the past 40 years. 

(A) If |si| < I is fixed, then, as x —>■ — 00 , 

w(x|s) = (—x)“3 1 /—2/3cos + /3ln ^8 (—x)®^ + + O ^(—, (1-18) 

with 

1 TT 

/3 = — In (1 - |sip) , (/) = --- argr(i/3) - argsi, 

and r(z) is the Euler gamma function. 

Remark 1.7. Expansion (1.18) appeared first in [1, 2] and was partially proven by Hastings, McLeod and 
Clarkson in [21, 10] using G elf and-Levitan type integral equations. In the late 80’s, Its, Kapaev, Suleimanov 
and Kitaev [23, 30, 27], using the “isomondromy method”, obtained the full leading term. Another derivation 
of (1.18) was given by Deift and Zhou [13] in the mid 90’s based on a direct asymptotic analysis of the 
RHP 1.4 with the help of the Deift-Zhou nonlinear steepest descent method [12]. We refer the reader to the 
monograph [19] for a detailed exposition of the methods used in the rigorous derivation of (1.18). 
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(B) If |si| = 1, then with si = ie, e = sgn(3si) G {±1} (compare here (1.11)), as a; —>■ —oo, 

M(a:|s) = (x"®) . 


(1.19) 


Remark 1.8. The leading order expansion in (1.19) is originally due to Hastings and McLeod [21] but was 
later on also derived by Deift and Zhou [12], who extended (1.19) to a full asymptotic series in which all 
coefficients can be computed recursively. 


(C) If |si| > I is fixed, then, as x —)■ —oo 
u(x|s) = 


with 


sin ^|(-x)2 + / 31 n ^8(-x)2^ + + O (^{-x) 2^ 

d= “ 1) ’ = -argr Q +i^^ -argsi 


+ 0((-x)-i) 


( 1 . 20 ) 


where the error terms are uniform outside some neighborhoods of the singularities of the trigono¬ 
metric function appearing in the denominator of the leading order. 

Remark 1.9. The singular asymptotic formula (1.20) was first obtained by Kapaev in [26] through the 
isomonodromy method. In [ 8 ], the authors rederived (1.20) based on nonlinear steepest descent techniques 
applied to the RHP l.f. 

In this paper we derive new asymptotic expansions of m(x|s) as x —> —00 with the help of a nonlinear 
steepest descent analysis applied to RHP 1.4. These expansions explain how the leading oscillatory power¬ 
like decay in case (A) is transformed to the leading power-like growth in case (B) as |si| f 1- We shall refer 
to the transition between regimes (A) and (B) as regular transition. In addition, we will also explain the 
transformation of the leading singular oscillatory power-like growth in case (C) to the leading power-like 
growth asymptotics in (B) as |si| f 1. It is natural to refer to the transition between (B) and (C) as singular 
transition. The results are as follows. 

1.1. Statement of results for regular transition. Let 

>r = ^ G ( 0 , 00 ); t = (—x)5 > 0 , 0 <u = — ln|l—|si[^|= — ln(l—|si|^) = —2^13, 0 < |si| < 1 ( 1 . 21 ) 
and k G (0,1) be implicitly determined from the transcendental equation 

2 k^ 


3 


l-bk2 


E' - 


1 -bk2 


K' 


( 1 . 22 ) 


in which K and E are standard complete elliptic integrals 


K = K{k) = [ 
Jo 


d/i 


A" = A'(k'); £ = E(k) = j( J 


E’ = E(k’) 

1 — /i^ 


/o v'(l-M^)(l-kV)^ 

with modulus k G (0,1) and complementary modulus k' = Vl — k^. It is shown in Proposition 3.2 that for 
any xr G (0, f'x/S), equation (1.22) determines k = k(xr) uniquely. In addition, set 




1 -bk2 




K 


T = r(>f) = 2 i—, 


(1.23) 


and let 


6 * 3 ( 0 , g) = 1-I-2 ^ q™ cos{2TTmz); 62{z,q) = 2 "^^ cos ((2m-b l) 7 rz), zGC 

m—1 m—0 

denote the Jacobi theta functions with nome q = The theta functions in turn determine the Jacobi 
elliptic function 


cd 20 AT 


l-k\ 1 -k 


1 + k / ’ 1 -bk 


6 * 3 ( 0 ,?) 02 iz,q) 

6 * 2 ( 0 ,?) 6 * 3 ( 0 , g)’ 


G C\ 


1 T 

- -b — mod Z -b rZ 


(1.24) 
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and we list properties of the elliptic integrals and theta functions which are relevant to our analysis in 
Appendices A and B. The main result concerning the regular transition asymptotics of M(a:|s) as x —> —oo 
is formulated in the following Theorem. 


Theorem 1.10. For any fixed S G (0, ^\/2), fi G (0, oo), there exist positive constants to = toiS),ti = 
= vi{fi) and cq = co((5),ci = ci(/i) such that 


i{x\s) = -e\ -- 


with 


and 


— , = cd ( 2(—x) ‘^Vix) K ( - — r I , - — r + Jl (x, s), 

2 y/\ + k’^{x) V Vl + kM/ 1+kW' 

\Ji{x,s)\ < cot~"^ Vt>to, 0 < X < t , 


2 2 

|-^i(a;, s)| <Vt>ti, -V2t - fi < v <-V2t. 

' ' in t 3 3 


(1.25) 

(1.26) 
(1.27) 


Theorem 1.10 describes the transition of the leading order asymptotics between case (A) for fixed |si| < 1 
given in (1.18) and case (B) for |si| = 1 in (1.19). Indeed, using the limiting behavior of V{x),t{x) and 
k{x) in the Jacobi elliptic function (1.25) as x I 0 and x f f (compare Corollary 3.3 below), we obtain 
directly 


Corollary 1.11. As x ^ —oo, 

u{x\s) = (—x)“3 COS ^8(—x)5^ + + O ^(—x)“to^ ^ (1-28) 

uniformly for 0 < x < . In addition, as x ^ —oo, |si| f 1? 

u(x|s) = + ^ 

uniformly for | ^/2 — \ < x < ^ \f 2 . 


1 

ln(—x) 


(1.29) 



Figure 2. Depiction of transition asymptotics. The region captured by Theorems 1.10 
and 1.14 is coined Boutroux region since the asymptotics are described in terms of Jacobi 
elliptic functions. Above and slightly below the separating line x = we observe to 
leading order Hastings-McLeod asymptotics, compare Theorems 1.12 and 1.16 below. In 
between we find Stokes lines, see Theorem 1.13. 

Our second result addresses the asymptotic regime x —?► —oo, |si| f 1 for x > fv^— “y with fixed /2 G M. 
In this case the leading behavior is already of type (1.19) and therefore connects the Boutroux behavior 
(1.25), (1.27) completely to the Hastings-McLeod asymptotics (1-19) for |si| = 1. 
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Theorem 1.12. Given /2 S ffi, there exist positive constants tg = tQ{f 2 ),vo = po(/ 2 ) and c = c(/ 2 ) such 
that 


i(x|s) = -e\ -- 1 - 


1 


Stt {tV2)2 


Y + J2(x,s) 


<J = - V^ — t = { — x) 2 

O 


(1.30) 


with 


yt>to, v>vo, v>-V2t-f2. 


\j 2 {x, s)| < ci ^ 

Theorems 1.10 and 1.12 describe the behavior of u(a::|s) as a: —>■ —oo for the values of x in 


0 , -v^ - 5 


U - y,oo^ C (0,oo). 


The corresponding regions in the double scaling diagram, i.e. Boutroux region for Theorem 1.10 and Hastings- 
McLeod region for Theorem 1.12, are shown in Figure 2. Another region, the Stokes region, corresponds to 


X S 


(^v/2 - (5, ^^2 - y) , (5 G ^0, , /i e (0,oo). 


(1.31) 


We will not provide a full asymptotic description of u{x\s) for the latter values of x as the nonlinear steepest 
descent analysis becomes increasingly difficult. Instead we focus on the scale 

t>to,v>vo- x>‘^V2-fJ^, /a e K, 

which eventually motivates the term Stokes region: For /a G (—oo, |) the leading order behavior of M(a;|s) 
turns out to be unchanged from (1.30), however once we cross the Stokes lines 


Sk ■■ v=-V2 t -—Int, k G Z>o, (1-32) 

once we go deeper into (1.31), additional {k + 1) terms will contribute to the leading behavior. This is in 
sharp contrast to Theorems 1.10 and 1.12 where the fixed choices of /i G (0, oo) and /2 G K had no effect 
on the leading orders. In the present paper, we shall prove the following estimation. 

Theorem 1.13. Given /a G (—oo, |) there exists positive constants to = to{f 3 ),vo = voifs) and c = c(/a) 
such that 


m(x|s) = 


• 1 + ^ ' 


ep 


+ Ja(x,s), 27rV2(2y2t)i 


(T = - V^ — X 
o 


(1.33) 


with 


|<^3(a:,s)| < ct A’il yt>to, v>Vo, V > -V2t — fslnt. 

(j 

This rigorously clarifies the appearance of the hrst Stokes line 5i: Fix /a G (—oo, |) in (1.33), then 

and subsequently by geometric progression, as a; —)■ —oo, |si| f 1: 

M(a;|s) = ^ ^ uniformly for —/a-y ■ 

1.2. Statement of results for singular transition. In this case we also work with the double scaling 
parameter 

XT = ^ G (0, oo); t=(—a;)5>0, 0 < = — In |l — |sip| = — In (|sip — l) = —27r/3, 1 < |si| < 

but notice that opposed to (1.21), we have now 

0 < — In (|sip — l) = —27r/3 = v ^ —27r/3 = — In (l — |sip). 

Furthermore we choose an upper constraint 1 < |si| < which is important to our analysis below, but 
at the same time is also consistent with the anticipated transition asymptotics for |si| ), 1.^ The module 


^For fixed |si| S [\/2, oo), the leading order behavior of ^(xls) as x —> —oo is already known through (1.20). 
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k £ (0,1) is again determined via (1.22) and we require V = V{>c) and r = t{>c) as in (1.23). However, 
instead of the Jacobi elliptic function cd(z) in(1.24), we need 


dc 2zK 


l-k\ 1 -k 


1 +k/’ 1 +k 


^ 2 ( 0 , g) 03{z,q) 

SsiO.q) d2{z,q)' 


2 ;GC\|^ mod Z + rZI . (1-34) 


The analogue of Theorem 1.10 for the singular transition is contained in the following Theorem. 

Theorem 1.14. For any fixed 6 £ (0, there exist positive constants to = to{S) and cq = co(d) such 

that 


X 1 + k(>!r) 


u(a;|s) = —e.i/ —— - ,- 

2 yi + k 2 (x) 


dc ( 2(—a;) = K 


1 — k(>!r) \ 1 — k(>f) 


1 + k(>£’) y 1 + k(>!:) 


+ Jiix, s) 


with 


I J4(a;, s)| < Cot 15 Vt>to, 0 < u < t , 

provided (a;, si) is uniformly bounded away from the discrete set 


Zn. = 


|(a:,Si): 2(—a;) 2 !/(>!:) = n £ Z\{0}| . 


(1.35) 

(1.36) 

(1.37) 


On the other hand, for any fixed fi £ (0,oo), there exist constants ti = ti(/i),ni = Vi{fi) and ci = Ci(/i) 
such that (1.35) holds true with 

2 2 

I J4(x, s)| < Vt > ti, V >Vi, -V 2 t — fl < V < -V 2 t 

' ' In t 3 3 

and no further constraint placed on {x, si). 

Here, Theorem 1.14 describes the transition to leading order between case (C) for fixed |si| >1 in (1.20) 
and case (B) with |si| = 1 in (1.19). Indeed, through the limiting behavior of V{x) and t{>c) as k f 0 and 
Kf in (1.35) we obtain 


Corollary 1.15. As x —>■ — 00 , 
M(a;|s) = 


sin ^|(—a :)2 + /31n (8(—a;) 2 ) + + O (^(—a:) 


+ (!l((-x)-w) (1.38) 


uniformly for 0 < k <t 5 and away from the zeros of the trigonometric function appearing in the denomi¬ 
nator of the leading order. In addition, as x ^ — 00 , |si| f 1, 


u(a;|s) = -ej-- +0 


1 


ln(—a;) 


(1.39) 


uniformly for | — j < >c < ^\f2. 

Our last result is the direct analogue of Theorem 1.12 to the singular transition. 


Theorem 1.16. Given /2 £ K, there exist positive constants to = to(/ 2 ),iio = co(/ 2 ) and c = c(/ 2 ) such 
that 

with 

2 

I ^ 5 ( 3 :, s)| < ct“^ Vt>to, V>Vo, V>-'/2t—f2- 

O 

The regions in the (t, z;)-plane captured by Theorems 1.14 and 1.16 are again shown in Figure 2 and we use 
the same terminology, Boutroux and Hastings-McLeod. We do not address the Stokes region for the singular 
case, a full analysis for all Stokes lines (1.32) in both regular and singular transition will be postponed to a 
forthcoming publication. 
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1.3. Further applications and outline of paper. The transition asymptotics are derived through an 
application of the nonlinear steepest descent method [12, 14] to the initial RHP 1.4. This analysis depends 
crucially on the values of the double scaling parameter 


V 



(0,oo); 


t = (—x) 2 >0, 0<'!; = — ln|l — jsip 


-In (1 - jsip), 0 < jsij < 1 
-ln(|sip-l), l<|si|<v^. 


For instance, the proof of Theorem 1.10 consists of several steps. 

(1) First, the master RHP 1.4 is analyzed with the constraint < v < — S),r] G (0,1) in 

place. All details are worked out in Sections 3 and 4, which include a g-function transformation, the 
construction of local model problems via Jacobi theta and Airy functions and iterative solutions of 
singular integral equations. This part resembles in some of its aspects the approach carried out in 

[7]. 

(2) Second, we analyze RHP 1.4 subject to the constraint t > > 0,t > tg with k G Z> 3 . The 

techniques used now are very different from (1), in fact they are close to the ones we would use in 
the asymptotic analysis of RHP 1.4 for fixed |si | < 1, i.e. no need for a g-function transformation but 
different model problems involving parabolic cylinder functions. We work out the necessary details 
in Section 5 and combining error estimations from (1) and (2), we obtain (1.25) with estimation 
(1.26). In addition, expansion (1.28) follows as well. 

(3) Third, estimations (1.27) and (1.29) are derived by applying modular transformations to the Jacobi 
theta functions used in (1) and introducing a new model function in a vicinity of the origin. The 
details are presented in Section 6. 

After that, we address the regime t > tg, v > vg, v > 1-^/2^ —/2 in Theorem 1.12. Here the nonlinear steepest 
techniques are again different from the ones used in the derivation of Theorem 1.10, we use a new gr-function 
and different model problems, compare Section 7. It is worth mentioning that these techniques resemble the 
ones used in the derivation of (1-19) with jsij = 1. Our final result for the regular transition is contained 
in Theorem 1.13 which we derive in Section 8. The appearance of Stokes lines is already observed earlier in 
Section 6, see Remark 6.4, however we give a quantitative analysis of this phenomenon only subject to the 
constraint 

i > to, V > Vo, /a e K. 

In [7], Theorem 1.12, the authors also discovered a Stokes region while analyzing P(s,j) as s —^ 00,7 f 1, 
recall Remark 1.3. Starting from the known large s asymptotics of the eigenvalues of ATgin all Stokes lines 
were classified by an application of Lidskii’s Theorem to the Fredholm determinant P(s,'y). We will now 
argue that we can reverse this approach and derive asymptotics of the eigenvalues {Aj(a:)}^g of ATai via 
Theorem 1.13 and identity (1.4). First, we state a slight improvement of Corollary 1.2. 


Corollary 1.17. As x —>■ —00 and 7 t 1? with cg as in (1.8) and a in (1.33), 


det (/ - 7A:Ai) 
uniformly for 


{ic,oo) 


= exp 






3 ■ t 

where b = 6(7) is positive, bounded in 7, such that 

b<a)^cg, 7tl- 


This result is a direct consequence of Theorem 1.13 and we state its proof in Section 8.2. Compared to 
Corollary 1.2, an additional contribution to the leading order appears in (1.40) and the extra factor is in 
general not close to unity. In other words, also the Fredholm determinant displays an asymptotic Stokes 
phenomenon, and we have identified the first Stokes line, 

2^1 

Si : V = -v2t- Inf. 

3 2 
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Now back to the spectrum of i^Ai, we let {Aj(a;)}^o denote the eigenvalues of the trace class operator 
i^Ai : oo); dA) O- These have been, for instance, analyzed in [31], were it was proven that the 

spectrum of Kai is simple, we have 1 > Ao(x) > Ai(ai) > ... and for fixed j G Z>o, 

\j{x) ~ 1 , a; —)■ — 00 . (1-41) 

In fact, Tracy and Widom in loc. cit. also gave an asymptotic formula. 


l — = 3^‘(l + o(l)), t = (—a;)2 —>■+00, 


J - 


(1.42) 


but the derivation is not fully rigorous. We will now outline a proof for the stated behavior, using Corollary 
1.17. Choose t > to,v > vq such that 


/a e ^- 00 , 0 ; z; = -ln(l- 7 ). 


(1.43) 


Apply Lidskii’s Theorem, 
det(/ - 'jKAi) 


det(/ - KaO 


L^(x,oci) 


= det {I + e-^KAiil - Km)-^) 


(XjOCi) 


nb 


j=0 


„-v Aj(a;) 


1 - Xj{x) 


= 1 + e 


-V Mx) 


1 - Ao(a;) 


^ det (/ + e-"ifi(/-ifi)-i) 


L^{x,oo) 


(1.44) 


with Ki = Kai • Pai and is the projection onto the eigenspace of Kai with eigenvalues {Aj}^;^. This 
allows us to compare (1.7), (1.40) to (1.44), i.e. as x —>■ —oo, |si| f 1 subject to (1.43), 

Ao(x) 


1 + e 


1 - Ao(x) 
though v 


det{l + e-^Ki{I-Ki)-^) = /"Mil + (l + o(t 

V Co 24 w+^co j ^ ^ 


-min{§,f-/ 3 } 


We multiply though with the second summand in the first factor in the right hand side of the last estimation, 

2 ty+tUo(x) ^ _ %) 

Co 


^ ^ ^.det{l + e-^K,{I-K,)-^) = MU (l + 2t (1.45) 

— Xo{x)) J Co V / 


and here, both factors in the left hand side of (1-45) are positive. But for t > to,v > vq such that 


/gelo 


7 


^ 2 /- . Int 

<+ d < a = -V2- K < — 


all summands in the right hand side of (1-45) are bounded. Thus, by positivity, all summands in the left 
hand side of (1.45) have to be bounded. So, using also (1.41), 

det{/ + e-”A'.(/-A'.)-‘) = 0(1), 

ea_ Xq [ x )) 

det {I + e-'^Ki{I - Ki)-^) = 0(1). (1.46) 

Since det(J + e~'’Ki{I — Ki)~^) > 1, we deduce from the first estimation, 

- = 0 ( 1 ), 

ei^Il-Ao(x)) 

which is consistent with (1.42) for j = 0. In fact, (1.45) allows us to deduce the inequality, 

1 - Ao(x) > + o(l)), (1-47) 

but in order to achieve equality our analysis in the Stokes region has to be extended; so far (1.46) only gives 


elV 2 t ^^ — Ai(x)) 


= o{i), 
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which is not sufficient yet to deduce equality in (1.47). Summarizing, provided Theorem 1.13 is extended 
to the full Stokes region and simultaneously also Corollary 1.17, the results on the transition asymptotics of 
M(a;|s) would enable us to derive (1.42) rigorously. 

The derivation of Theorem 1.14 follows largely its regular counterpart. 

(4) The master RHP 1.4 is analyzed asymptotically subject to the constraint < v < t(|-\/2 — <5), rj € 
(0,1) in Sections 9 and 10. Opposed to the regular case (1) we require a different outer parametrix 
which has the singular structure and corresponding exceptional set (1.37) encoded. 

(5) After that, we address the scale t > 1 ;^+^ > 0,t > to following [8]. The details are summarized in 
Section 11. 

(6) Again, by modular transformations, plus an additional argument related to the singular structure, 
we complete the proof of Theorem 1.14 in Section 12 and in addition obtain Corollary 1.15. 

Section 13 concludes the manuscript by deriving Theorem 1.16 and we list a few identities for complete 
elliptic integrals and Jacobi theta functions in Appendices A and B. 


2. Preliminary steps for regular and singular transition analysis 

We start with the scaling transformation Ar(A) = Y (^Xy/—x) , A S C\ (Ji and are lead from the initial 
RHP 1.4 to a RHP for the function A (A). This problem is formulated on the same jump contour IJ® Pj, with 
jump 

A+(A) = A S P* 

and jump exponent 

^^(A) = i (^A^ - a) , t = {-x)Y 

Next we deform the original jump contours, to obtain a RHP for a function Z{X) as shown in Figure 3 - at 
this stage the endpoints A = ±A* are yet to be determined. 

Riemann-Hilbert Problem 2.1. The function Z{X) has the following analytical properties: 

• Z(X) is analytic for X € C\([—X*, X*] U Uf jk) 

• We have the following jumps, compare Figure 3, 

Z+(A) = Z_(A)e-‘’’(^)'^^5fee‘’’(^)'^^ A e 7 fc, A: = 1,..., 6; 

Z+{X) = A e [-A*,0]; 

Z+{X) = A e [0,A*], 

where in fact 

= a2S3S4S5a2, ^2 = 

• As A —)■ 00 , the function Z(X) is normalized as 

Z{X) = I + o {x-^). 

Next, referring to the domains Tlj,j = 1,..., 4 shown in Figure 4, we set 

'g-tii(A)cT3 5--igt,?(A)<T3^ AeHi 

g-U?(A)<T3 5--lgU5(A)<T3^ AGH2 

T(A) = Z(A) <1 A G H 3 (2.2) 

AGH 4 

/, else, 

and are lead to a RHP for T{X) with jumps (compare Figure 4 for the jump contours) 
r+(A) =^_(A)e-*’’(^)‘"^5'fee‘’’(^)'"^ A G 7fe, fc = l,3,4,6 

T+(A) = ^_(A)e-*’’(^)'"^5'4e*’’(^)'"^ A G 72; T+(A) = A G 75 


0 -i 
i 0 


( 2 . 1 ) 
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Figure 3. Deformation of 
original jump contours, 
y(A) i-A Z{X). The three 
dots represent the points 
A = -A*,0,A*. 


72 



Figure 4. First opening of lens, 
Z{X) hA T(A). 


and 


T+(A) = A G 72 ; 


T+(A) = ^_(A)e-*’’(^)'^^CT 252 ^ 4 ■V 2 e‘’’(^)'^^ A G 75 
r+(A) = T_(A)e-‘’’(^)‘"='o•2S'4S'35'4CT2e‘’’(^^'"^ 


for A G [—A*, 0] and A G [0, A*] in the last two cases. Notice that 




1 ~ S1S3 Si + Si(l — S1S3) 
-S3 1 - S1S3 


0'25'4S'3>S'4(T2 — 


1 - S1S3 


-S 3 


Si + Si(l - S1S3) l-SiS3y’ 


(2.3) 


The presence of the factors 1 — S 1 S 3 = 1 — |sip on the diagonal motivates the use of a 5 -function in order 
to address the transition asymptotics x —)■ — 00 , |si| —>■ 1. Subsequently we shall first work out the necessary 
details in case of the regular transition. 

3. Regular transition analysis for xt g [^, |a/2 - >5] with 0 < <5 < ^^2 fixed 
We introduce the double scaling parameter 


In (1 - Isil 


^ t 


-3 -- -^G(0,oo); t=(-a;) = , /3 = ^ \n {l - \sif) (3.1) 

(—X)2 t Z-K 

and keep x G [<5, § 7/2 — (5] with 0 < (5 < 57/2 fixed throughout. Based on the cyclic constraints (1.11) the 
parameter x allows us to derive expansions of the Stokes multipliers in terms of x which are used later on. 

3.1. Preliminary expansions. First observe from (1.11) and (3.1) the exact identity 3?(si) = 

Second, by definition (3.1), 

|si| = 1 —-I-O (e“^^*) , t —>■ 00; arg(si) = earccos ( | , arccos : 1,1] —^ [0, tt], 

2 V Dll / 

and we obtain for Si = 

Proposition 3.1. As t ^ 00, |si| f 1 with = ie, 

Si = ie ^1 — -(1 -I- ies 2 )e -1- O (e ^ sg = —ie ~ ies 2 )e + O [e ^ 

and 

Si + Si(l — S 1 S 3 ) = ie ^1 -f -(1 — ies 2 )e + O [e , 

uniformly for x G [ 5 , §7/2 — < 5 ] with 0 < 6 < 57/2 fixed. 


(3.2) 


(3.3) 
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Let US now move ahead and introduce one of the key ingredients in the nonlinear steepest descent analysis, 
the ^-function. 


3.2. Introduction of g-function. Set J = (—1, — k) U (l,k) C M where k € (0,1) is determined implicitly 
in (1.22), 

[ \/(l-M^)(M^-k2)dAi. (3.4) 


3 V 1 + k2 


2 k^ 

E' - - -Twif' 


l + k 2 ’ 


1 + k 2 


The following Proposition addresses the bijective correspondence between x € (0, fv^) and k S (0,1). 
Proposition 3.2. The modulus k is uniquely determined via (3.4) in case x G (0, Iv^)- Moreover, as 

X 


1 1 n 2xr 29 /XT 

k = l-2W- +-— (- 

V TT TT 




and, as X \ with a = a{x) = — x I 0, 


k = , 

Proof. Consider the function 
where 


I Incrl 


ln|lncr| 2 + 7ln2 

1 + ^ + —n-+ ^ 

2 In O’ 4 In cr 


F(k, >f) = >r-I(k), ke[0,1], 


In I In a\ 
In^ cr 


(3.5) 


(3.6) 


I(k) = 


l+k 2 


A/(l-M^)(M^-k 2 )d^. 


With the help of standard expansions for the complete elliptic integrals as k 0 and k f 1 (compare Appendix 
A), we get 


I(k) = ^72 


o o on 

1 - -k^llnkl - -(1 + 41n2)k2 + —k^] lnk| + O (k^) 

24 16 ^ ' 


k I 0, 


I(k) = ^(1-k)^ 


l + (l-k) + ^(l-k )2 + 0((l-k)3) 


, kfl. 


(3.7) 

(3.8) 


Thus 


limF(k, >c) = >c-< 0, limF(k, >c) = >c > 0, 

k^O 3 ktl 


but F(k, xr) is real analytic in a neighborhood of (k G (0,1], x: G [0,oo)) with the first partial derivatives 
equal to 


Fk(k, xr) = k 


1 + k 2 


l+k^ A 


\/(l - - k2)d/x + 


cm 


d/r 


> 0 , 


F^(k, xr) = 1. 

Hence the implicit function theorem guarantees existence of a unique real analytic solution k = k(xr) of the 
equation F(k, xr) = 0 near the point (k, xr). We use (3.8) and (3.7) to derive (3.5) and (3.6). □ 

Besides the expansions (3.5) and (3.6) for k = k(xr) itself, we will later on also require expansions of the 
frequency V = F(xr) and module r = r(xr) introduced in (1.23). These follow directly from (3.5) and (3.6) 
and are summarized in the Corollary below. 


Corollary 3.3. As x f 0, 
2 


= -to “ + *" ° = -(1" (i^) 

Secondly, as a = — x fO, 


X \ 17i xr 


87 r TT 


V{x) = -- 


Incrl 


ln|lncr| 2 + 71n2 „/ln^|lncr| 

1 + H- 

In cr 2 In cr V In cr 


+ 0{xi). (3.9) 


(3.10) 
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and in the same limit also 


^_^ |lna| 

r(>f) 27ri 

We now define the g-function, 


Inllncrl 7 In 2 Inllncrl 

1 - 1 -^ ^- -+^ 

In tr 2 In cr In a 


In^ a 


g{z)=4i - m^))^ dfj., zeC\[-M,M] 

Jm 


where 


0 < TO = 


< 1 <M= ^ ^ 


< 


1 


(3.11) 


(3.12) 


(3.13) 


V^v/TTl? 2 y2vT+T? V2 

and the contour of integration is chosen in the simply connected domain CP^\[—M, M]. Moreover, we fix 
-TT < arg ((^2 - — TO^)) < TT such that \/{gP — — m?) > 0 for /r > M, 

and thus, g{z) is single-valued and analytic in CP^\[—M, M]. As z tends to infinity 

'l-k 2 ^ " 




+ 0 (z ^) , z —)• oo, 


where 


x/(Al2-M2)(^2-TO2)-/r2+_ 


d/ 1 . 


(3.14) 


e = -d(M) + 41 [ 

JA 

Further steps in the analysis require certain analytical properties of g{z). 

Proposition 3.4. For z G M introduce 

= '^{ 9 +{z) + g-{z)), and n(z) = g+(z) - 5 _(z), with 5 ±(z) = hm 5 f(z ± k). 

The functions fl{z) and n(z) are real-valued on the real line, in fact 

n(z) = —8 f J (/l2 — M2) (/l2 — 77 j 2 ^ z€{M,+oo), ^(z) = 0, zG{m,M), 

J M * 

n(z) = —8 f ij (M 2 — /i 2 ) (r7j2 — /t 2 ) d/i, zG(—to,to), 

J Z 

n(z) = —8 f \J (m 2 — fi'z'j (m 2 — /i 2 ) d/t = 27rP (>f), z G 

J —m 

n(z) = —8 f (m 2 — g?') (to 2 — g?) dfj, + 8 j ij (/i 2 — M 2 ) (/i 2 — m 2 ) d/i, zG(— oo,—M). 

^ —m J z 


Moreover 


n(z) = 0, z G (—oo, —M) U (M,-l-oo), n(z) = 8 / ^(M2 — g?') — m^) d/i, zG{m,M), 

J Z 

n(z) = 8 f \J (M 2 — /i 2 ) (^/i 2 — d/i = x:, zG(—to, to), 

J m 

n(z) = XT— 8 / (M2 — /i2) (/t2 — 77 i 2 ^ 2;G( —M, —to). 

J z 

3.3. The g-function transformation. We first go back to the RHP for T(A) as defined in (2.2) with jump 
contour shown in Figure 4. Now fix the endpoint 

A* = ^ =M 

72 Vl + k 2 

and employ the following transformation, 

,5(A) = A G C\St, Sr = M] U [J U 7 U 7 U 72 U 75 ^ , (3.15) 
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with g = g{z) as in (3.12) and £ in (3.14). Since all jumps in the T-RHP display the structure 

r+(A) = T_(A)e■‘’^(^)‘"=^GT(A)e‘’’(^^‘"^ A S C\Sr, 
we are lead to the following RHP 


Riemann-Hilbert Problem 3.5. Determine the 2x2 matrix valued function S(X) = S{X;x,s) such that 

• <5'(A) is analytic for X € C\Ss where the jump contour Eg is identical to the contour 'St shown in 
Figure 4- 

• The function S{X) satisfies the jump condition 

S+{X) = S-{X) e-‘9-(^)'^3G.J,(A)e‘9+(^)'"^ A G Eg = St- 


=Gs(A) 


As X 


S{X) = I + 0{X-^) . 

Recall at this point (2.3): for A G {m,M) we have thus 

with K from (3.1). But (compare Proposition 3.4) 

XT—n(A)=8 f ^ (M^ — — TO^) d /1 > 0, X£{m,M) 

J m 

and n(A) > 0 in the right slit. Hence, using also the expansions (3.2),(3.3), 




-1 


J, X —>■ — CX), |si I t 1 • ^ 




0<5<-^/2 (3.16) 

O 


uniformly in A chosen from any compact subset of the right slit In the left slit (—M, —m) a similar 

situation occurs, 

Q (\\-q /XX /e"‘(^"n(A)) _ g^gg))gita(A)\ M \ 

‘^+(A)-^-(A)(^_53e-i‘f2(A) (1 _ s^53)e-*n(A) j. AG(-M,-m). 

>^-n(A) = 8^ ^ — m^) d^ > 0, A G (—M, —to) 


But here 


and also 


n(A)=8y > 0, Ag(-M,-to). 

Hence combined with (3.2),(3.3), 


GsW 


0 

_g —i^e—itn(A) 


gif e+itn(A) 

0 


-1 


J, a; —>■ —oo, |si| t 1 : G 


(5,-y2-(5 


, 0 < (5 < 

O 


(3.17) 

uniformly in A chosen from any compact subset of the left slit (—M, —to). Next, we consider the jump 
contours which extend to infinity, i.e. 

Ssoo = 72 u 75 U 72 U 75 U IJ 7fc. 

Along these contours the jumps in the S'-RHP are given by 

S+{X) = S'_(A)e-‘9(^)'"='G^(A)e*s(^)'^^ A G Es^,^, 


Hence, by triangularity and the sign chart of 5 ft( 5 (A)) (compare Figure 5), the jumps on the infinite contours 
approach the identity matrix exponentially fast in the limit x —> —oo, |si| f 1 with k G [d, Ia/S — d], d > 0, 
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provided we stay away from the endpoints A = ±to, ±M and the origin A = 0. In the remaining gap 
{—m,m) C K two cases need to be distinguished: First for A G (0,m), 


s+W = 


and secondly for A G (—to,0), 


S+{X) = 


f 1 —\ 

V(si + si(l - SiS3))e"**^(^) (1 - sisa)^/ 

/ 1 0 \ (1 

S_{X)Sl,WSu,W, AG(0,m) 


(_s 3 ei*f 2 (A) 

•5-(A) 

5_(A)5l,(A)5c,,(A), 


(si + Si(l - siS 3 ))e‘‘^('^)A 
(1-S1S3)^ J 

0\ /l (si + Si(l - SiS 3 ))e“^(^) 

l) 1,0 1 

A G (—TO, 0). 


(3.18) 


(3.19) 


In (3.18) and (3.19) all off-diagonal entries are fast oscillating as x —>■ —oo, see Proposition 3.4. We now 
transform this behavior on (—to, to) to exponential decay with the help of contour deformations {opening of 
lens) tailored to the factorizations written in (3.18) and (3.19). The key to this explicit transformation is 
the following Proposition. 


Proposition 3.6. Introduee for z G {—m, to) the functions 

Hi{z)=in{z), H 2 {z) =-ifl{z), 

with n = n(z) as in Proposition 3.4- Then Hi{z) admits local analytical continuation into a neighborhood 
of the gap (—to, to) into the upper half-plane such that 

fi{Hi{z)) < 0, > 0, iRz G (—TO,to). 

Similarly, the function H 2 {z) admits local analytical continuation into a neighborhood of the gap (—to, to) 
into the lower half-plane such that 


fi{H2{z)) < 0, Xsz < 0, ftz G (—TO,to) 

Proof. We follow the standard line of argument. First 

in(z) = —8i f \J (M^ — p?') {m? — p'^') dpt, z G {—m,m) 

J Z 

and thus 



II 

0 

II 

1 

00 

K5 


II 

0 

II 

1 

00 

K5 


which implies the claim via the Cauchy-Riemann equations. 


□ 


3.4. Opening of lens. Suppose = 1)2 denotes the lens shaped regions shown in Figure 5. With the 
help of the analytical continuations of Hi{z) and H 2 {z) as discussed in Proposition 3.6, we set 


L(A) = 


S{X)Sfjl{X), 

S{X)SlA^), 

S{X)Sf}^{X), 

5(A)5i,(A), 




A G 

XGCf 
X G C2 

A G £2 

otherwise 


(3.20) 


so that L{X) solves the following RHP 
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3i(<?(A)) > 0 


72 



Figure 5. Opening of lens, the jump contours for L(X) consists of the solid black lines. In 
addition, along the dotted lines 3?(g(A)) = 0. 


Riemann-Hilbert Problem 3.7. Determine the 2x2 piecewise analytic function L{X) such that 

• L{X) is analytic for A G C\ {[—M, —m] U [to, M] U U 7 ^ U 

• We have the following jump behavior, with orientation as indieated in Figure 5, 


• A —>■ 00 , 


L+(A)=L_(A) 


Gs(A), 

X G i-M, 

^£7, (A), 

XG^ft 

Sl,W, 

X G 7 f 

Su,W, 

A e 7 ^ 


A e 7 ^ 

L{X) = 

i + o{x- 


) U (to, M) U 


After employing the explicit transformation (3.21), it is now time to focus on the local model problems 
near the points A = ±to, ±M and A = 0 as well as on the slit segment J = (—M, —to) U (to, M). 


3.5. The outer parametrix. The RHP associated to the outer parametrix is motivated by (3.16),(3.17) 
and thus consists in 


Riemann-Hilbert Problem 3.8. Find a 2 x 2 matrix-valued piecewise analytic function N{X) = N{X; t, e) 
such that 

• JV(X) is analytic for X G C\ J 

• Along the branch cuts, with orientation as in (5), 

/ Q if£+itn(A)\ 

N+{X) = iV-(A) 0 Ag(-M,-to) (3.21) 

N+{X) = N_{X)(^J., (3-22) 

• N{X) is square integrable on J 
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• A —>■ OO, 

N{X) = I + O {X-^) 

Observe that the function 

7V(A) = A e C\J 

solves a similar RHP as the one posed for N{X), but with jumps on J given by 

7V+(A) = A^-(A)(^_^_;;n(A) q Ae(-M,-m), 

iV+(A) = N_{X)(^^^ Ae(m,M). 

This is precisely the same jump behavior which appeared in [7], Section 3.1. The solution method is therefore 
analogous, we only summarize the relevant steps and refer to [7] for further details. 

The solution to RHP 3.8 is derived in terms of Jacobi theta functions defined on the elliptic curve 
P = {{z,w) : = p(z)}, p{z) = {z"^ - m^) - M^) 

of genus one. We view T as two sheeted covering of the Riemann sphere, glued together in the standard way. 
For definitness, let yjp{z) ~ as z —)■ oo on the first sheet and y^p(z) ~ —z^ in the similar limit on the 
second sheet. 





Hi 


di 


Figure 6 . Standard homology basis for F 


We fix a homology basis for F as indicated in Figure 6, and let 


c dz / X i 


•y/ (M^ — p?){p?' — m?) 
be the unique holomorphic one form on F with normalization 

w = 1, 


-1 


i M 
2 ~K' 


and H-period (compare (1.23)) 


T = ® LO = 2c 
Ibi 


dz 


= 2i- 


K 


— IT > 0. 


OO 


(3.23) 


/_„ - z^){M^ - z^) K'^ 

Also, define the Abel (type) map 

u : CP^\[—M, M\ —>■ C, z H-)■ m(z) = f uj = u{oo) - \- O {z~^) , z ■ 

Jm ^ 

and collect the following properties 

Proposition 3.9. The Abelian integral u{z) is single-valued and analytic for z G CP^\[— M, M], we have in 
addition 


u+(z) — U-{z) = 


0, z S (—OO, —M) U (M, oo) 
— 1, z S (—m, to) 


, u+(z)+M_(z) = 


0, z £ (to, M) 

T, z £ ( — M, —to). 


and 


u{oo) = - (in CP^). 


(3.24) 
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Remark 3.10. Equation (3.24) is a simple identity between complete elliptic integrals, 

K 


i K 

u{oo) = - 


2K'^ 


Our construction requires furthermore the functions 
’ {z + m){z — M) 


uj{z) = 


{z + M){z — m) 


> (t>{z) = \ (w(z) + {uj{z)) , (^(z) = ^ (uj(z) - (a;(z)) , 


defined and analytic for z G CF^\[—M,M] such that a;(z) > 0 for z > M. To complete the derivation, we 
use the Jacobi theta function 


and define 


0(z|t) = 03 {z\t) = ^ exp [iirk'^T + 27rifcz] , z G C, 

fcGZ 

^ / 0iuiz)+tV±d) 0i-uiz)+tV±d) \ ^ ^ 

^ ’ d{u{z)±d) ’ 0 {-u{z)±d) ^ ^ ^ ’ 


where (see (1-23)) 

V{k) = V = -/ \/ (M^ — ^ 2 ) {jji2 _ ^ 2 ^ = 0 ( 2 )^ 2 ; G (—M, —to); 

^ J — m 

With this, a solution to RHP 3.8 is given in the next Proposition. 

Proposition 3.11 (cf. [7], Section 3.1). The function 

N(X) = pif-aiR ( N^-^\X)UX)\ -G 

0{tV) i-iv(-)(A)<^(A) Nt\X)cf{X) 




(3.25) 


(3.26) 


is single-valued and analytic in C\J. Its jumps are stated in (3.21) and (3.22), furthermore, as X ^ 00 , 

. e'jtv) 


1 


NiX) = I+- 


—c 


e{tv) 


X \ e{ 0 ) e{u{oo)+tV-d) M-m ^-i^e 


9{u{oo) — d) 


2i 


e( 0 ) e(u{oo)-tV-d) M-m 

6{tV) 9{u{oo) — d) 2i 

9'{tV) 

^ e{tv) 


o{\ 


- 2 ^ 


Our next move focuses on the construction of model functions near the branch points and the origin. 


3.6. Parametrix near the origin. Near the origin, we first observe that 

2g{X) = X — iTt{X), A G 72 ; 2g{X) =—x — iTt{X), A G 75 . 

Thus the jumps on the quasi-vertical segments 72 U 75 displayed in Figure 5 approach the identity matrix 
exponentially fast along the entire segment, more precisely 

g-t9(A)<T3 5.^5.-lgtg(A)a3 ^ J + + S 2 )e-‘^e‘‘^(^V+ , A G 72 

= /-(si + S 2 )e-‘"e-™(^V_, A G 75 - 

Using also (3.2) and (3.3) we shall consider the following RHP. 

Riemann-Hilbert Problem 3.12. Find a 2 x 2 piecewise analytic function H{X) such that 

• H{X) is analytic for A G (C n D{0,r))\{'jf U 7 ^) where 0(0, r) = {A G C : |A| < r} and we fix 
0 < r < f. 

• The boundary values are related via the equations (compare Figure 5 for orientation) 

R+(A) = H_{X)[^^ j, AGO(0,r)n(7i+U72+); (3.27) 

R+(A) = R_(A) (^.^^_L(a) J), a G 0(0, r) 0(71-0 72”) (3-28) 

• i7(A) is bounded at the origin A = 0. 
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• As X ^ —oo, |si| t 1 with >c G [<5, — <5],S (0, |\/2) fixed, we have 


uniformly for X G dD{0, r). 


H{X) = {l + o{l))N{X) 


The solution to this problem is very elementary: assemble 


= 






pK<^3 


1 iee**^(°) 


,0 


1 


1 


O' 


_igg-ita(o) I 


argCG 

argCG (f,f) 


argC G 


with r2(0) = ttV (see Proposition 3.4) and observe that H°{() solves a “bare” RHP 

_ 3 ir ' 


H°{C) is analytic for C G C\{C : argC = ±f ,±^} 
Along the four rays oriented from zero to infinity, 


1 

0 1 


HliC) = HliC) 


, argC = 


1 -iee**^(°) 

0 1 


, argC = 


i?^(C) = ffl(C) 


Stt 

argC = -j] 

Stt 

, argC = --. 


• H°{() is bounded as C —0. 

• As C —>■ oo, we have 

H-°(C) = (/ + C>(C”°°))e‘‘i'"3 (3.29) 

uniformly in a full neighborhood of C = oo. 

Referring to the locally analytic change of variables, 

C(A) = ^(H(A)-H(0)) = 4tMmA(l + C>(A2)), A S D(0,r), 0 < r < ^ (3.30) 

the origin parametrix is then given by 

Pr(A) = Af(A)i^°(C(A))e-'‘i(^i'"^ XGD{0,r), (3.31) 

with N{X) as in (3.26). Since N{X) is analytic in the disk D{0,r), compare Proposition 3.11, we check 
directly that the jumps of H{X) are indeed as in (3.27) and (3.28) with the orientation of the contours near 
the origin like in Figure 5. Note also that the four rays in the bare RHP can always be locally deformed 
to match the contours Furthermore with (3.29) and (3.30) we obtain from (3.31) the desired matching 
between the model functions: as t = {—x)i -G +oo, |si| f 1 such that >f e [d, f-s/S — d] is fixed, 

H{X) = {l + 0{t-^))N{X) (3.32) 

uniformly for 0 < ri < \z\ < r 2 < y- This completes the construction of the origin parametrix. 

3.7. Parametrices near the inner branch points. Our construction for the model function near the 
right inner branch point A = m is motivated by the local expansions 

iH(A) = co(A — to)2 + CA ^(A — m)2^ , A S yj*" H iA(m, r) 

—iH(A) = —co(A — m)2 + CA ^(A — to)2 ^ , A G y(" n iA(77i,r) 

n(A) = XT — Co(A — m)5 + CA ^(A — m)5^ , X G {m, M) n D{m,r) 


where 


Co = ^-\/2 to(M2 — TO^) > 0, 

O 
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Figure 7. Jump Figure 8 . Jumps of the model 

contour for the bare function C/(A) near X = m. 

parametrix {Q. 


we fix 0 < r < min{^, — m)} and the function (A — m) 5 is defined for A G C\[to, oo) with the branch 
fixed by the requirement arg (A — m) = tt for X < m. This motivates the use of Airy functions and our 
construction follows from now on [7], Section 3.2. with minor modifications. We define 


A«"(C) = Ao(C)<^ 


^f -f^ 
^1 -f^ 

where Ao(C) denotes the unimodular entire function 


argC G (0, 
argCG (f,f) 


J N, argC G (f^,27r) 


(3.33) 


Ao(C)=iv^e 


0 



Ai (e ‘" 3 ’' Ai(C) 

'^Ai'(e-‘^C) Ai'(C)^ 


e‘ 3 ‘^ 3 ^ CeC 


(3.34) 


which is constructed with the help of the Airy function w = Ai(z), the unique solution to the boundary 
value problem 


w = zw: 


74 9 A 

Ai(z) = —= 

Zy TT 




("■’) 


00 , — TT < argz < TT. 


Through the standard properties of w = Ai(z) (cf. [29]), the model function {() in (3.33) has jumps on 
the contour shown in Figure 7, more precisely 

• A^^(C) is analytic for C G C\{argC = 0, 

• The boundary values are related via the equations 

A 


A""(C) = A""(C) 


A""(C) = A""(C) 


0 1 

1 0 

-1 1 


27r 

«gC = 


, argC = 


Jtt 


A«^(C) = A«^(C)(^_\ J), argC = 0. 


As ^ > 00 , we have in a full neighborhood of infinity, 

1 f-1 6i 




With the (locally conformal) change of variables 


48C^ 


6i 1 


+ 0 






(3.35) 


C(A) = 


6teA {{M^ - fi^) {m^ - dfi 


2 

{At\/2m{M‘^ — m?) ^ ^ (A —m), |A —m| < r, (3.36) 
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we define 

C/(A) = |A - to| < r 

which involves the (locally) analytic function 


(3.37) 


i3,,(A)=7V(A)e'5' 


1 1 


5(A)-3 C(A) 


X-M\^ 


, 5(A) = 


X — m\^ 
X-M 


1 . 


A — m _ 

Remark 3.13. Analyticity of Br^{X) near X = m follows from the simple observation that for X G (rn,m+r), 




0 


—1 1 

1 1 


(<5-(A)) 


‘"3g-ifc73 


C(A) 


X — m 


= {Br,W)_, 

thus Bri{X) can only have an isolated singularity at X = m. But this singularity is at worst of square root 
type, hence has to he removable. 

As we are allowed to locally deform the jump contours in (3.33) and match those in the initial L-RHP, 
the local parametrix U{X) has jumps in D{m,r) as shown in Figure 8 . Moreover, with the help of (3.35), 


U{X) = Ar(A)e‘A"‘"3 
= I + N{X) 


48(5 V^i 1 
1 f-l - 6 e 
48(1 \ 6 e 1 


+ 0 


(c-i) 


o[c"^)\{NiX)) ' 


N{X) 


so that (as t —>■ +oo, |si| 11 with k g [5, f's/S — 5] fixed) 

U{X) = {l + 0{t-^))N{X), 

uniformly for 0 < ri < |A — to| < r 2 < min{™, |(M — to)}. 

The parametrix near the remaining inner branch point A = —to is constructed similarly: introduce 

(I, arg( e (-7r,-|) 

arg( e (-f,f) 


(3.38) 

(3.39) 


O' 


A^^{C) = AoiO{ 


— g-i7r(l-7) I j ’ 

" 1 0 \ 


(3.40) 


, arg(e(f, 7 r) 


where 


with 


^ ^-g-i7r(l-7) ij \0 1 

Of pm 

7 = 1H-/ ^(M^ — ^2) _ 2tV(x), 

^ J — m 


AoiC) = iv^e‘5(i-T')‘-3g-if<T3 


X gif'^3g-if(l-7)^3 


0 A /e'’^Ai'(e‘’^() —e* 3 Ai'(e* 5 (^ 

0 e~'^J \ —Ai (e‘'^() Ai(e‘5(^ 


(gC, 

and collect the following properties 

• is analytic for ( G C\{arg( = — |, |, 7 r} 

• The jumps are as follows, compare Figure 9 for orientation. 


(3.41) 


^^"(0 = A(!"(() 
I^^(() = A^^iC) 
A^^iO = A^^iO 


1 0 

_g-i7r(l-7) 2 ^ 

I gi7r(l-7) 

0 1 

X gi7i-(l-7) 

_g-i^(l-7) Q 


, arg( = -- 


, arg( = 


, arg( = 7 r 
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argC = I 



Figure 9. Jump 
contour for the bare 
parametrix 


-itn(X) 


Figure 10. Jumps 
of the model function 
F(A) near A = —m. 


1 

0 1 



1 0 

i£e-“U) 1 


As C —>■ oo, valid in a full neighborhood of infinity, 

1 


RH(^\ _ -i7r(l-7) 




-1 -ie'’"(i-T') 


48C5 




O 


(ri) 


63 


jiC^cra 


(3.42) 


With the change of variables 


C(A) = 


6t {{M^ - ^i^) {m^ - d^l 


2 

^ (A + m), |A + TO|<r, 


we move ahead and define the parametrix near A = — to, 

F(A) = |A + to| < r. 


(3.43) 


where 


B,,(A) = lV(A)e‘?' 


_igi7i-(l-7) jgi7i-(l-7) 


1 


A + M\ ^ 


with J(A) = ^ 1, A ^ cx). 

Remark 3.14. Again, the multiplier i?£j(A) is analytic at X = —to, since 
{Bi,iX)),=N.{X) 


0 gif e+27i-ity\ /_igi’r(l-7) jgi’r(l-7) 


/+ —\'V ^ _ ^—i^e—2'7TitV 


e 4' 


>-(A))- 


/ \ /Vf \ “ 4‘^3 

X |^C(A)y 3 p^j ={Be,{X))_, Xe{-m-r,-m). 


Moreover, we check that F(A) (after employing a local contour deformation) has jumps inside D{—m,r) 
as shown in Figure 10 and through (3.42), 


F(A) = N{X)e'^‘ 


/+fg. 


48C3 


i) + 0 (r®) 


e 3 


iV(A). 


Thus for t —)• + 00 , |si| t 1 such that >c € [J, |\/2 — J], J S (0, jV^), 

ViX) = {l + 0{t-^))N{X) 

uniformly for 0<ri < |A + to| < r 2 < min{™, |(M — to)}. This completes the construction of the 
parametrices at the inner branch points A = ±to. 


(3.44) 
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3.8. Parametrices near the outer branch points. Once more Airy functions are used in the construc¬ 
tions, first near \ = M: Define the bare parametrix, 


A^^iO = MO< 



with 


^o(C)=o'2Ao(C) (T2=V7’‘e‘5 


7 = 1 


-1 1 


0 

o-if 


argC G (-TT,-^) 
argC G (-^,-f) 
argC G (-f, f) 

argC G (f,7r) 


Ai(e‘3^) Ai (e^’^C) 

e'5Ai'(e'5C) e*'^Ai'(e'^^C) 


(3.45) 




This model function has the properties listed below 

• A«^(C) is analytic for C G C\{argC = f,7r} 

• We have the following jump behavior on the contour shown in Figure 11, 

^^"(0 = ;), «g(=-^ 

j). »g< = » 


argC = f 



Figure 11. Jump 
contour for the bare 
parametrix A^^(C). 



Figure 12. Jumps of the model 
function P{X) near X = M. 


As C —> 00 , 


I' .i J 

This time, we use the change of variables, 

6 t [ - M^) - m^))^ dn 

J M 


48(5 


-1 -6i 

-6i 1 


O 


(o 




(3.46) 


C(A) = 


{X-M), \X-M\<r 
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and the local parametrix is defined as 

P{X) = |A - M| < r. 


(3.47) 


Here, 


Br,{X) = N{X)e'^' 


1 -1 

—i —i 


C(A) 


A — m A '* 


<5(A) = 


X-M\^ 
X — m 


1, A —>■ c» 


X-M, 

is once more analytic near A = M. The relevant jump properties of P{X) are displayed in Figure 12 and 
with (3.46), 


P(A) = iV(A)e‘T"'"^ 
= I + N{X) 


48C 


48C 

i f-1 6 e 
- 6 e 


;")+0(rS)|(Ar(A))-' 


N{X). 


Hence for t —>■ +oo, |si| t 1 such that x e [d, — <5],<5 G (0, ^v^), 

P{X) = {l + 0{t-^))N{X) 
uniformly for 0 < ri < |A — M\ < r 2 < \{M — m). 

We are left with the parametrix near the left most branch point A = —M: assemble 

/1 

(o 1 )• 

/ 1 o\ A 

L-i7r(l-7) ijlo 


(3.48) 


H«"(C) = Ho(C)<^ 


1 


argC e ( 0 , f) 
argCe (f,f) 
argCG(f,f) 


, argC e (f^,27r) 


(3.49) 


where 

^(0 = 




, CgC, 


■e‘f 0 \ /Ai'(C) -e '3"Ai'(e H)\ 

0 e 'ey Ai(e“'^^) j 

and, as before, 7 = 1 — 2tV{x). We have 

• A^^iC) is analytic for C G C\{argC = 0, f, 

• Fixing orientations as shown in Figure 13 below, the relevant boundary values are related as follows 

,0 1 


H«"(C) = A«"(C) 


TT AtT 

: arg C = -, — 


H«"(C) = A«"(C)( , argC = 


271 


H"" (C) = A"" (C) 


yir(l-7) 


, arg C = 0 


As C —t 00 , 


Igif 


-i 1 
i 1 


„-if (1-7 )ct3 


— 1777 


1/1 6 ie ' 

48^ -1 


+0 (c-s) 


g-§C^-3 


(3.50) 
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Figure 13. Jump 
contour for the bare 
parametrix {(,). 


Figure 14. Jumps of the model 
function Q{X) near X = —M. 


and with the change of variables 


C(A) = 


- M'^) - m^))^ 

J-M 


the required local parametrix is given by 

Q(A) = \X + M\ < r, 


(4ty2M(M2 - m2) ) ' (A + M), \X + M\<r, 

(3.51) 


where 


Bi-AX) = Af(A)e'5'^'^3gif (i-7)^3 


1 -1 

—i —i 


g-if (1-7)<T3A(A)-'"3 [ 


X + m\ ^ 
X + M 


and (5(A) = { ^ —>■ 1 as A —>■ oo. Since Ufa (A) is analytic at A = —M, the model function Q{X) has the 

jumps shown in Figure 14 and with (3.50), 

1 f 1 6ie-‘’^T'^ 


Q(A) = 7V(A)e‘4-3 
= I + N{X) 


48(f VSie'"'" -1 


O 


(<-*) 


48(5 


1™ ®'‘ir>°(f-") 


— 6 ee 


iV(A). 


This implies for t —>■ +oo, |si| f 1 such that k G [(5, f — <5],(5 S (0, |\/2), 

Q{X) = {l + 0{t-^))NiX) 
uniformly for 0 < ri < |A + M\ < r 2 < — m). 


(3.52) 


3.9. Final transformation - ratio problem. In this final step, we use the explicit model functions H{X) 
in (3.31), U[X) in (3.37), F(A) in (3.43), P{X) in (3.47), Q{X) in (3.51) and N{X) in (3.26). We transform 
the previous RHP 3.7 for L{X) by setting 


R(A) =L(A) < 


{H{X)) \ 
(UW) \ 
{v{x))-\ 

(PW) \ 

(QW) \ 

{N{X)) \ 


|A| < r 
\X — m\ < f 
|A + m| < r 
\X-M\<f 
\X + M\<f 

|A| > r, |A =F w| > f, \X =F M| > f 


(3.53) 


where 0<r< ™,0<r< min{^, |(M — m)} and 0 < f < 5 (M — m) remain fixed. Definition (3.53) 
implies that the RHP for L{X) shown in Figure 5 is transformed to the RHP for R{X) posed on the contour 
Yifi shown in Figure 15. 
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72 



Figure 15 . The jump contour Ejj of the ratio function i?(A) defined in ( 3 . 53 ). 


• i?(A) is analytic for A S C\Ei^ where = [—M, —m] U [m, M] U yj*" U 7^ U Cq U U U Ci^ U 

U 

• For the oriented contour shown in Figure 15 , the jumps read as follows. First inside the five 
circles Cj, 

i?+(A) = i?_(A)Gfl(A;t, |si|), Gfl(A)S_(A) = 5 -(A)G(A) 

where 


E{X)=H{X), |A|<r; 



|A — m| < f 
|A + m| < f ’ 


S(A) 


P(A), \X-M\<f 
Q(X), \X + M\<f 


and (inside the far most left circle Ci^), 


(\ (ie —Si)e 2ts(A)\ 

lo 1 


G(A) = 


y(s3 + ie)e2*sW ij’ 

I -iesi(l + 

I _jgg-itn(A)g-i(x+n(A)) 


_jg gitn(A)g-t(>.-n(A))' 
ies3 


A S (72 U 74) n D{—M, r) 
A G 73 n D{—M, f) 

X e i-M, -M + f) 


followed by (inside the left circle G^J, 


G(A) = 


1 (si + Si(l - S1S3) - 

^ 1 o\ 

-(s3+ ie)e-'*^(^) li’ 




—iesi(l + e (si + si(l — S1S3) — ie)e‘"“""'‘'e 

— 2 t>c 




le S 3 + e 


A G 7^ n D{—m, r) 
A G 7^ n D{—m, f ) 
X G (—m — f, —m) 
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as well as for the corresponding circles on the right, 

r/l + 


,0 


1 


C(A) = 


1 , 0 , 
^(si + Si(l - S1S3) - ie)e Ij 

ie S3 


_i£g-*(^-n(A)) 


_igg-t(>r+n(A)) _igs^(i + e"^*) + , 


A S 7^'' n D{m, f) 
A G 7f n D(m, f) 
A G (m, TO + f) 


and 


C(A) = 


1 0 \ 

^(si — ly 

''1 —(s3 + ie)e“^*®(^A 

_igg-d^-n(A)) 


A G (71 U75) n D{M,r) 
A G 76 n D{M, f) 


le S 3 


igg-t(><+n(A)) 


Inside the remaining circle centered at the origin, 

' (l (si + S2)e-"'‘e“(^A 

0 1 y ’ 

1 o\ 

-(si+ S2)e-"^‘e-“(^) ly’ 

1 -(s3 + ie)e‘*^(^)\ 

0 1 y ’ 

1 (si + sie-"'* - ie)e'*^(^A 

0 1 j 

1 i)\ 

-(S 3 + ie)e-“(^) ly’ 

1 0^ 

I (si+sie-^*-ie)e-“(^) 1 . 


C(A) = 


A G 72 n D{ 0 ,r) 
A G 75 n D{ 0 ,r) 
A G j+nD{ 0 ,r) 
A G 7^ n 11 ( 0 , r) 
A G 72" n l?(0,r) 
, A G 7f n £>( 0 ,r) 


Recalling our previous expansions for si = S3 in ( 3 . 2 ),( 3 . 3 ), we see that for A G Sr inside the circles 
with radii r, f and f, we have 


C(A) = / + C> (e-^‘) , t^+(X),|si|tl: e 


( 5 ,- 72-5 


which is uniform with respect to the chosen radii. But since .^-(A) is bounded on the same contours, 
we obtain the estimation 


\\Gr{-; t,\si\) - I\\L2nL°°{^anDj) < die di > 0 


( 3 . 54 ) 


where Dj denotes any of the five open disks centered at A = 0 or A = ±to, ±M. Secondly we 
estimate the jumps on the remaining finite branches of Sr: for the horizontal line segments, first for 
A G (—M + f, —TO — f), followed then by A G (to + f, M — f). 


R+(A) = R_(A)iV_(A) ( '' 


_• itn(A) -t(j<-n(A)) 


R+(A) =R_(A)iV_(A) 


lee 

ies 3 


le S 3 


(fV_(A)) \ 


_• -t(>f-n(A)) 


_ig g— *(^+n(A)) _• 


■si(l + e-^*) 


( 7 V_(A)) \ 
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Notice that for A G {m + f,M — r), 

pX p\ 

k-u{x) = 8 — m^) d/i > 8-\/2m(M + m) / — m) d/i 

J m J m 

> c-\/k(l — k) f2, c>0 universal, 


with a similar estimate also holding for A G (—M + r, —m — f). As the outer parametrix N_(A) 
remains bounded on the line segments for fixed radii, we have with universal dj > 0 (using again 
(3.2),(3.3)), 


|si|) - /|U.nL=e((-M+.-.-m-f)U(m+GM-.-)) < 


(3.55) 


Remark 3 . 15 . Although we have chosen fixed radii in (3.53) we shall always indicate the dependency 
of the estimations of Gfi{-]t, |si|) — I on the latter. This will be useful later on when we partially 
drop the constraint on >c G [J, — 5]. 


Next, on the lens boundaries, for j = 1, 2, 


R+{X) = R_{X)NiX)Su,{X){NiX)) \ XGj+n{\XTm\>n, 
R+{X) = R.{X)N{X)SL,W{N{X)y\ XGj-n{\XTm\>f}, 


and therefore, with fixed radii which ensures the boundedness of N{X), 


l|Gfi(’ |si|) L^nL°°{-yfn{\\^m\>f}) 




(3.56) 


The pieces of Hr now left are all infinite branches and circle boundaries, more precisely 
R+(X) = R_{X)Q{X){N{X)y\ X G Ce,- R+{X) = R_{X)V{X){NiX)y\ X G C,, 


i?+(A) = i?_(A)17(A)(iV(A)) \AGai; -R+(A) = R-(A)P(A)(Af(A)) \ A G a. 

and 

R+{X) = R.{X)H{X){N{X)y\ AgCo. 

We conclude with (3.39),(3.44),(3.48),(3.52) and (3.32), (fixed radii so that N{X) is bounded) 

I|Gr(- |si|) - d|U2nL~(uc,) < y (k(l - k))~" max|f"5, Vkf-ij (3.57) 

Also, for the infinite branches, 

i?+(A) = P_(A)Af(A)e"‘9(^)'"3GT(A)e*®(^)'"3(Af(A))"\ 

where Gt{X) depends only on the Stokes multipliers and Gt(A) are piecewise constant triangular 
matrices. Due to this triangularity and the sign chart of 3?((/(A)), compare Figure 5, we have that 

||G,^(-;t,|si|)-/||L=nL~(infi„ite) d, > 0. (3.58) 


• As A —>■ 00 , 

R{X) = I + O {X-^) . 

The latter RHP for R{X) can be solved iteratively as t ^ + 00 , |si| f 1 provided >c G [d, ^^2 — <5]. 
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3.10. Iterative solution. We collect estimations (3.54)-(3.58), so that 


IIGrI'; |si|) ~ ^l|L2nL“(EH) < t —>■+00, |si| t 1 : G 




0 < (5 < ^V2. (3.59) 

O 


Here we used that for the values of G [5, ‘^y/2 — i5] we only work with hxed radii in (3.53). Since the RHP 
at hand is equivalent to the singular integral equation 

R_{X) = I+^ f r_(w)^Gr{w)-I)^^, AgSr, 

2711 ' w- X 

estimation (3.59) guarantees (cf. [12]) the unique existence of an asymptotic solution of the RHP in L'^{T,r), 
this solution satisfies moreover 

kil) - I\\L^i^R) < (3.60) 

3.11. Transition asymptotics. We now extract the asymptotics through (1.16), i.e. we apply 

u{x\s) = 2 lim A(y(A;x, s)) 

and for this have to trace back the sequence of transformations 

y(A) hA a:(a) ha z{x) ^ t(a) ^a s{x) ^ l{x) ^ r{x). 

First, we get that 


m(x|s) = 2 lim A(P(A)),„ = 2\/—x lim [zX(z)l = 2y/—x lim [ 2 :T(z)l 


= 2 ^/—x lim 

z—>-oo 

= 2y/—x lim 


2—>-C50 


12 

= 2 \/—x lim 

J 12 2 —>-oo 


J 12 






12 


But since 


we have 


Also, as z —)■ 00 , 


g[z) - •d{z) = £ + 


l-k 2 
8 z Vl + k 2 


+ O (z ^) , z —>■ 00 , 


g-t(g(z)-i5(z))cr3 _ g-ttcrs / j_ 


it /l-k^ 


8z V 1 + k2 


a, + 0 {z-^)]. 


R{z) = I+ 7 ^ [ R-iw){GRiw)-l)^^=I+^ [ R.{w){GR{w)-l)dw + 0{z-^) 

JSb W — Z ^11^ JEr 


and 


iV(z) = /+ —+0(z-2), z- 


with (compare Proposition 3.11), 


iAT\ ^(0) 0{u{oo) - tv - d) M - m & 3 (Olr) 0 2 {tV\T) ■^^y M - m 

[Iyi)l 2 = ~7:7TTTT-;:7—;—^-^r;-^ ^ ^ —r—rC -;r;-e 


0{tV) 9(u{oo) — d) 2i 

where we used (3.24) and (3.25), i.e. 


02(O|T)d3(tP|T) 2i 


w(oo) -d = 


and the well known transformation identities between Jacobi theta functions, compare Appendix B. So far, 
we have thus the exact identity. 


i(x) = V—3 


„2U 


— e{M — m) 


6 * 3 ( 0 |f) 02 {tV\T) 
02{o\t) e^itviT) 


+ S{k) 


(3.61) 


with 


£{>c) — ~ I (j^-{w){Gr{w) — l)^ dw = O (t , t —>■ +00, |si| t 1 : G 


(5,-y2-(5 
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In order to further simplify expansion (3.61), we note that 

Proposition 3.16. 

Proof. Let us treat 

r 

f{M) =e + ijy = -d(M) + 4i / 

2 7a 

— 2i / m = m{M) = 

J — m ' ^ 

as a function of one variable M G [ 5 : H is not hard to verify that 


A/(/r 2 - M2) (^2 _ ^2) - ^2 ^ _ 


d^ 


—/(M) = 2iM(m2 - m2) 


d/i 


d/r 


IM 


(^2 _ - m2) 7-m -^Z(M2 - ^ 2)(^2 _ ^ 2 ) 


But from the identity m(c») = compare (3.24), we see that the last difference of integrals vanishes, hence 
/(M) is in fact M independent. Letting M ), d, we have 


/(M) = ^+o(l)-^ = o(l), Mi^ 
and thus the stated identity follows. 

Summarizing our previous simplifications and noting that £ G iK, we have 

u(a; s = -ev^(M-TO) ’ + Ji{x,s), t = -x)^ 

6'2(0|r) 03 (tV\T) 

and 


□ 


(3.62) 


Proposition 3.17. For any given 6 G (0, ^v^) there exist positive constants to = to{S),vo = vo{S) and 
c = c(S) such that 

|7i(x,s)| < yt>to{6), v>vo{6), tS < V < t ^-V2 — 6^ . 

The leading term in (3.62), although currently expressed in terms of ratios of Jacobi theta functions, can 
be rewritten with the help of the Jacobi elliptic function (1.24) using Appendix B, 


9j{z\t) = 9j{z,q), 7 = 2,3; 
and the identity (see (3.13)) 


ITTT 

q = e = exp 


^ K 


.r 1 1-k 

M — m = 


Hence, 


m(x|s) = -eJ-- 


X 1 — k 

2 ym? 


72 Vl + k2 ■ 

l-k\ 1-k 


cd 2{-x)^VK 


l+kj’1 +k 


+ 7i(x, s) 


(3.63) 


and the error term Ji(x, s) is estimated in Proposition 3.17. We have already control over the error in the 
domain 


I (t, u) : t > to, V > Vq : tS < v < t |, 


but in order to complete the proof of Theorem 1.10 this domain will be extended in the following sections. 
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4. Extension at the lower end for regular transition 
A ssume that both t = (—x)® > to and v = — ln(l — |sip) > Vq are sufficiently large such that 

0 < < -v^ — (5, with 0 < 77 < 1 fixed. (4-1) 

3 

From (3.5) and (3.13), we see that 

as t > to for some c > 0. Hence, in the definition of the ratio problem (3.53) for i?(A), we cannot keep the 
radii f and f fixed since M — m ]. 0. This issue can be resolved by scaling both radii with t, in fact we shall 
choose in (3.53) 

f = Cit“3, f = C 2 t~^, 0 <Ci,C 2 : Cl + C 2 < c. (4.2) 

In order to estimate the jumps in the ratio problem with the latter choice in place, we need to derive 
estimations for the outer parametrix N{X) and A^-(A) on the jump contours. First with Corollary 3.3, 
expansion (3.9), 

03 (z|r(>f)) = 1 + C>(>f), >f|0 

uniformly for z chosen from compact subsets of C. Next, with the contracting choice of radii, 

uj(z) = 0{1), = 0{1), 0 < — (5 

o 

uniformly for z G Cg^ U Cg^ U U Cr2- The same estimations for w^^(z) are also valid on the infinite 
branches as well as on all remaining finite branches outside the five disks Dj centered at A = 0 and 
A = ±m, ±M. Together 

iV_(A) = 0(l) = (iV_(A))-\ >ri0 

uniformly on the circle boundaries and all other finite or infinite branches outside IJ Dj. We can now simply 
go back to (3.54)-(3.57) and substitute (4.1),(4.2) and (3.5) into the estimations. We obtain immediately 

llGfl(-;t, |si|) -/||L 2 nL=e(s^) < t > to, v > vo : Q < < x <‘^V2 - 5, 77 G (0,1), (4.3) 

and thus, repeating the steps in the previous section (with adjusted error terms), 

»(xW = -^-fl ^ c<l (2(-.)iv' A- (1^) , + <2 ((-«)-'U.) , (4,4) 

uniformly as x ^ — 00 , |si| f 1 such that 0 < < x < — S for any 0 < 77 < |. In estimation (4.4), we 

obtained control over the error term through the a-priori knowledge 

p_(-;t,|si|)-/|U2ni=o(s„)<^ (4.5) 

and a direct estimation of 

= O ((-a;)-^+5'7^ . 

We can get a better error estimation by explicitly computing the first terms in 
= {Gniw) - I)^^dw+i (Gniw) - I)^^dw 

+ [iGR{w)-l){R-{w)-l)'j^^dw. 

Two of the four contour integrals along the circles Gg^UGr^ have been explicitly computed in [7], Section 4.2. 
It was shown in loc. cit that the aformentioned two contour integrals are in fact O even with contracting 
radii. The same result also applies to the contour integrals over Gg^ and and without reproducing the 
lengthy computations of [7], we simply conclude 

^ {{Gniw) - I) {R-{w) -l))^^dw + 0 (t-i) . 


(4.6) 
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Remark 4.1. Even without referring to the explicit computations in [7] we see direetly from, say (3.36) 
and{3.38), that the evaluation of the contour integral 


j> {Gr{w)-I)Aw 


by residue theorem leads to an asymptotie series in reciproeal nonnegative integer powers oft. 

For further improvement, we use that we have slightly better -estimations than (4.3), 

||Gfi(-;^, |si|) - I\\l^{^r) < 

t 4^/ 

SO that by general small norm theory of [12], 

|si|) - i\\l^{^r) < • 

t 4^ 

Back to (4.6), for any r/ € (0,1), 

£:(>f) = c> +o{t-^) 

and we have therefore derived (3.63) with an error term Ji(x,s) such that 

Proposition 4.2. For any given 5 € (0,^y/2),p G (0,1) there exist positive constants to = to{S,r]),vo = 
Vo{6, rf) and c = c{5, rj) such that 

|•^l(a^, s)| < yt>to, V > Vq, t^~^ < V < t — 6^ . 


Notice that for t > to such that 0 < x: < t 5 we obtain from (3.5) and (3.9), 

l-k\ 1-k 

TTk/ ’ TTk 


—e \/ — ? ^ cd ( 2(—x) K 


2 cm? 


(y—x) 4 \/—2/3 cos{'KtV (>r)) -b O (^—x) 10 ) 


with 


1 


= ^ln(l - |siP) ■ 


Also from (3.9), as 0 < x: < t 5 , 

TTtV{x) = — /31n ^8(—a;)5^ -I- /31n | ln(l — |sip)| — ,5(1 -|- ln27r) -b O ((—a;)”™) ■ 

But from Stirling’s approximation, 

argr(i/3) = argF = “/3(1 + ln27r) -b ,51n | ln(l - |sip)| -b ^ -b O ((xrt)“^) , 

since xrt = 0(ti) —)• -boo. Using in addition (3.2), 


fir 

arg(si) = y+ C)(i”~) 


expansion (4.9),(4.8) combined in (4.7) gives therefore 

"2 

7rtU(xr) = — 
where (compare (1-18)) 


(-a;)2 -/31n ^8 (-x) 2^ -b 0 - ^(1-b e)-b O (^(-x) 1 ^°) 


(4.7) 


(4.8) 


(4.9) 


(4.10) 


(j) = --- argr(i/3) - argsi. 

Substituting (4.10) back into the expansion of the Jacobi elliptic function, we use the addition theorem for 
cos z and the identity 

TT 

—ecos (2 -b — (1 -b e)) = cosz, z S C, 











TRANSITION ASYMPTOTICS FOR THE PAINLEVE II TRANSCENDENT 


35 


SO that finally 

(^) ’ r^) V^cos 0(-a;)i + /31n (8(-a:)i) + 

+ C>^( —x)“5^, X —>■ —OO, (4-11) 

uniformly for 0 < k <t~'s. This short computation shows that the Jacobi elliptic function leading term in 
(3.63) reproduces the known oscillatory leading behavior of (1.18) for the values of {t,v) such that t > to 
and 0 < >r < <“ 6. However in the same region of the double scaling diagram we do not yet have control over 
the error Ji(x,s), compare Proposition 4.2. This will be achieved below using different nonlinear steepest 
descent techniques applied to the RHP 1.4. 

5. Further extension at the lower end for regular transition 

Suppose throughout that t = (—x)5 > to is sufficiently large such that 

t > >0, t>to, k G Z>o. (5.1) 

With this constraint in place we employ a different approach to the nonlinear steepest descent analysis of 
RHP 1.4. In fact we now choose steps which are close in certain points to the ones used in the derivation of 
(1.18) for fixed |si| < 1, see [19], chapter 9, §4. 

Start from the Riemann-Hilbert problem for Z{X) with its jump contour shown in Figure 3 or Figure 16 
below and where we fix 



The jump matrix on the full segment [—A*, A*], compare (2.1), is factorized 

s;'sGs;' = .ASA,, = (‘" ^ J) e-- (J ‘) s 

and we observe that with (5.1), 

0 < < exp , t>to, t> > 0. (5.2) 


SR(t9(A)) > 0 



Figure 16. Opening of lens subject to the scale (5.1), the jump contours for <i)(A) are the 
solid black lines. In addition, along the dotted lines 3?('d(A)) = 0. 
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Next define (compare Figure 16), 


ci>(A) 


A € rif/ 

^ (A)^^, A s ri/, 

Z{X), otherwise 


and obtain the following RHP 

• $(A) is analytic for A S C \(72 U 75 U IJ® 7 ^) 

• The boundary values are related via the following jump conditions 


$+(A) 


XGjk, k = 

$+(A) 


A G 72 

$+(A) 


A G 72 

$+(A) 


A G 75 

$+(A) 


A G 75 

$+(A) 

= ^-{X)Sd, 

A G [-A*,A* 

The function <i>(A) is 

normalized. 



$(A) =/-hO(A-i) , A^ 

■ CX). 


Notice that the jumps on the infinite branches 7 ^ for A: = 1,3,4 ,6 are exponentially close to the unit matrix 
in the double scaling limit (5.1) as long as we stay away from A = ±A*. On the remaining infinite branches 
we have for A G dD{\*,r) U dD{—X*,r),r > 0 with our choice A* = ^, 

G^(A) + 

with a universal constant c > 0. Thus recalling (5.2), we could choose a contracting radius r = r{t), compare 
Section 5.3 below, and obtain again sufficiently fast decay for the jumps on the latter infinite branches. We 
will work out all necessary details after the following two subsections. 


5.1. The outer parametrix. The required outer parametrix is given by 

AgC\[-A*,A*]; 1 .= ^ (5.3) 

and appeared in [19], chapter 9, §4. Its relevant analytical properties are summarized below 

• 4'^(A) is analytic for A G C\[—A*, A*] 

• The boundary values are related via the equation 

(A) = ^'(?(A)e-^*'"X A G (-A*, A*) 

• 41 ^(A) is square integrable on the closed interval [—A*, A*] 

• As A 00 , we have 

4>^(A) =/ + ^ +C>(A"2) . 

A 


5 . 2 . The edge point parametrices. Near the end points A = ±A* the required parametrices are con¬ 
structed in terms of parabolic cylinder functions, see again [19], chapter 9, §4. We briefly summarize the 
results: 


Let D^{() denote the unique solution to Weber’s equation 


w" + 


1 C"\ . 

--- li; = U, w 


w{C) 


subject to the boundary condition 

i^.(C)=re-H= 



K^- 1 ) 

2C^ 


+ o(r") 


37r 371 

--<argC<-. 
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Observe that w = is entire in ^ G C. We define for C G C the Wronskian type matrix 

■(u+l) qN 

0 1, 


7 rn - 9-1-3 (D-.-iiK) D,{0\ 0 

ZM) - 2 I d’M) 


and assemble 


z""(C) = Zo(C)<^ 


argC e (-f,0) 

argC G (0, f) 
argC G (f ,7r) 


where 


and 


-ffn = 


1 0 
ho 1 


H, = 


1 hi 

0 1 


I, 

Ho, 

HoHi, 

H0H1H2, 

H 0 H 1 H 2 H 3 , argCG(f,f) 

i /„+2 n = 0,l,2 


argC G (tt, 

f Stt Ttt > 


ho — 


\\f^ 


hi = 


1 + hohi = e 


27rii/ 


r(i + p)’ r(-p) 

Using standard properties of parabolic cylinder functions, cf. [29] this setup leads to the following RHP for 
the bare parametrix 

• is analytic for C G C\{argC = 0, |,7r, 

• The jumps along the contours shown in Figure 17 are as follows 

Z«^(C) = Z^^iOHk, argC= ^ = 0-1>2,3 


^(C) = argC=^ 


7 B.H/ 


As C 00 , valid in a full neighborhood of infinity, 




1 1 
1 -1 




{y) 2 s 

{v + l)2s - {v)2s 


{-Y{-v-1)2s ) s!2- 


X g(V-(‘"+5)'”0-3 

where we used the full asymptotic series 


^.(C) 


_ cl 
e 4 




(-0 


2 s .-2s 


Stt 


Stt 


s=0 


s!2^ 


C ^ 00 , —- < argC < — 


combined with the identity 


^lC) = |^.(C)-^.+i(C), C,^GC. 

Next, we employ the locally conformal change of variables 

C(A) = 2v^(i?(A*)-i5(A))^ =2V^e'^(A- A*) (^1 + ^(A-A*) 
and define the edge point parametrix near A = A* as 


|A-A*| <r. 




where 


«W=(c(A)^)'''”(-h)""”e-»<-)-.2-P.(«'') J), |A-A 


-A < r 


is analytic at A = A*. After a local contour deformation we check that the jumps of ^'’'(A) near A = A* are 
identical to the ones in the initial <I>-RHP, i.e. (compare Figure 18) 
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Hi 



Figure 17. Jump 
contour for the bare 
parametrix {(,). 



Figure 18. Jump contour of the 
model function 'I'’’(A) near A = A*. 


5-;(A) = A G 71 ^; ^\{\) ^, A G 72 ^ 

«';(A) = A G 76^; ^'^(A) = A G 75 ^ 

= AG(A*-r,A*), 


where 'jjr = 7j H D{X*,r). Moreover, with the previous asymptotic expansion of Z^^{Q 

h, \ 


S 3 


A 

1 _ 


vL 


®’-(A) ~ {/)(A))"> f- 

L C(A)^ V 

00 

X E(-r((-u)2,n-(-U-l)2An) 


CJU y. 

/ + ^ C2s{l^) 


— 2s 


s=l 
— 2m 


s!2" 


S 3 


(/3(A)) "^vii^(A) 




m—1 


/•/'\A —2m+l ^ 

<(A) ^ ^ 


m! 2’" 


771=1 


f (u)2m 

0 \ C(A)-2™1 

l 0 

(-)'"(-u)2+ m!2™ J 


^'^(A), (5.4) 


as t —>■ +00 (and thus |C| —t 00 ). Here we used the abbreviation 

^(A)= (c(A)^^)%-‘^("*). 


Expansion (5.4) establishes the matching relation between the model functions 4'’’(A) and 4'^(A) valid as 
t —>■ 00 for 0 < ri < |A — A*I < r 2 < j. 

Near the other edge point A = —A* we can use the symmetry of the problem and define the required local 
parametrix as 

5-^(A) = a2^’’(-A)a2, |A + A*|<r (5.5) 

which replaces the matchup (5.4) by 


4'^(A) 


+ (7) S*"+ - 


C(-A) 


— 2m 




C(-A) 


— 2771+1 


771=1 


d2" 


m! 2' 

00 

'E 

771=1 


T.{/3i-x)y 


(-)”^(-u)2m 0 

0 (u)2m 


C(-A) 


— 2771 ■ 


i!2’^ 


^^(A). 


(5.6) 


Collecting the model functions 4''°(A), 4"’(A) and 4/^(A) we move on to the final transformation. 
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5.3. Ratio problem. Assemble 


x(A) 


r(vI/’'(A)) \ |A-A*|<r 

$(A) (vI/^(A)) \ |A + A1<r 

[(vI/^(A))-\ |A±A*|>r 


(5.7) 


with 0 < r < j which we will choose more specifically below. This leads to the RHP formulated on the 
contours as shown in Figure 19 




Figure 19. Jump contour for the ratio problem as introduced in (5.7). 

• x(A) is analytic for A S C\{Cr U Q U {infinite branches 7 ^}) 

• The following jumps are valid: 

X+i\)=X-(.m^WSj{'f^{X)r\ As 7,, J = 

along the infinite branches 7 ^ with piecewise constant matrices Sj which are given in the <i>-RHP. 
There are no jumps inside the circles and along the line segment [—A*, A*] since we constructed 
“exact” local parametrices 4'^(A), ^'’'(A) and Ji^(A). However on the circle boundaries, 

X+(A) = x_(A)vI/”(A)(iE-^(A))-\ Aea; x+(A) = x_(A)iI/^(A)(vI/^(A))-\ AeQ 

• As A —>■ c», we have that 

x(A) = / + 0(A-i). 

For the jumps on the infinite branches 7 ^-, j = 1,..., 8 we obtain immediately 

Gx(A) = J + 

with universal constants di,d 2 > 0. If we choose a (in general) contracting radius r such that 

1 1 L 

- > r = - > t fc+i > 0, 

4 V ~ 

then by (5.1) for t > > 0,t > to, 

di>ct — d2tr^ < dif^ — c? 2 t^ < —dsf^, d^ > 0, k > 3; 
hence all contributions from the infinite branches are approaching the identity matrix exponentially fast, i.e. 

l!G;,(-;t,|si|)-J|U 2 nL~(uy,) fc e Z>3 (5.8) 

for t > > 0,t > to- For the circle boundaries we use (5.4) and (5.6), for instance for A G Gr we use 

|C(A)|>4y2t™ as well as 


hi=0 ((>ct)5et^‘) 
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which follows from Stirling’s approximation. Thus for X G Cr 


|G'x(A;t, |si|) - l\<c- 




0 




y/t 


Kt -ip(r) 


+ £(A;t, |si|) 


where 


! \ ^ f /-tW^ 

^M = -~glC(A)3^ 




< — arctan(r) = —r il + O (r^)) = — (1 + o(l)). 

TT TT ' ' ' TT 


Bounds for the error term f(A;t, |si|) can be found for instance in [29]: there exists to = io{k) > 0 and 
constants dj = dj(to) such that 

fc -1 

|lf(A;t,|si|)|| Vt > > 0, t>to, keZ> 3 . 

For the latter estimation to hold for fc = 3, we have to adjust the radius r i—>■ - with a sufficiently large, but 
k,t, jsij independent, constant c > 0. Notice that with (5.5) completely similar estimations will also hold for 
G^{X;t, jsij) when A € Ce. Thus, summarizing 

;t,\si\) — I\\L2r,L'=°(CrUCe) ^ diot ^(fc+i) ^ Vt > > 0, t>tQ, fc S Z>3. (5.9) 

The latter estimation combined with (5.8), we obtain from general theory [12] the unique solvability of the 
ratio problem for the values of {t,v) satisfying (5.1). Moreover its unique solution satishes 


\\X-i-',t,\si\) - < diit =('=+!), Vt>u*’+^>0, t>to, keZ> 3 . 


(5.10) 


5.4. Derivation of leading order asymptotics subject to (5.1). We need to trace back the transfor¬ 
mations 

r(A) hA X{X) hA Z{X) hA $(A) hA x(A) 

in order to obtain with (1.16), 

M(a:ls) = 2 lim A(F(A)),„ = 2-\/^ lim [2 <i)(z)] = 2-\/^ lim [z4'^(z)x(z)] 

A—>-oo L ' ' z—>-oo 2 —>-oo 


= Zy—x 


ly—x 


va3 + 


271 


X-{w)(G^{w) - I) dw 


<f {G^iw)-l)„dw + o({-x)-^), fcGZ>4 

JCrUCi ^ ' 

where we already neglected all exponentially small contributions from the infinite branches, compare (5.8). 
Also the latter error estimation follows from (5.10) and the following computation: observe that by residue 
theorem and (5.4) as well as (5.6), 

£ P^{w)C{w)-^^+^dw 

= (-2Ti)e-2*’’(^*) +0 (f-™) ; 


= (-27ri)e^*’’^'^*) + O (t ^ . 


Recalling the definition of hi and applying standard symmetrization techniques, we obtain with fc = 4, 

u{x\s) = {—x)~i \/—2/3cos (8(—a;)5^ + £ + O (^(—, (5-11) 
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uniformly as a; —>■ —oo subject to the scale (5.1) with fc = 4. Estimation (5.11) appears in Corollary 1.11, 
expansion (1.28). Moreover, with (4.11), we have now shown that 

/ 1 ^ I ^ 1 —k , — k\ 1 — k\ ^ , . 

“(*!•) = TttP (2(-»^) w k (rqrt j • TTk j 

where 

Proposition 5.1. There exists positive constants to and c such that 

\Ji{x, s)\ < ct~'^, Vt>to) 0<v<ti 

Propositions 5.1 and 4.2 combined together (for 77 = | in the latter) form the content of Theorem 1.10, 
estimation (1.26). In order to complete the proof of Theorem 1.10 we have to derive the outstanding 
estimation (1.27) by extending the region 

{t,v) : t>to'. 0 < — (5 

at the upper end. This will be achieved in the next section. 



6. Extension at the upper end for regular transition 


In the analysis of the lower end extension presented in Section 4 we used the crucial fact that for >«■ 0, 

we have —ir —> + 00 , compare (3.9). Thus at the lower end, the oscillations in the Jacobi theta functions 
become strictly periodic and their amplitudes vanish exponentially fast. In turn we obtained 

' ^ \ l\ \ iC > ' „ \ I ' 


N_{\) = 0{1) = \ >ri0, XgJ:r\\Jd, 


which lead us to the central estimation (4.3). At the upper end the oscillations in the theta functions still 
vanish, however this time not through vanishing amplitudes but increasing periods: as f we get with 
(3.10), 


271 

I Incr 


+ =0(1), a=‘^V2 -KlO. 


This implies the necessity of using modular transformations for the Jacobi theta functions appearing in the 
construction of the outer parametrix (3.26), i.e. we use (cf. [29]) 


/ |^/\ ^ITTT Z 


= \j - 03 {zT'\T')e 


02(2|r) = \ - 94{zT'\T')e''' 


, z e C, r' =-. 

T 


Going back to (3.26) the transformation identities imply for A G C\ J, 


IV(A) =e'5' 


, ^(OkO „-2.u'tu.(oo).3 ( N[^^\x)^iX) iV'(+)(A)^(A)\ 2.ix'iu„(A).3e-il 

Oir'tVlT') \-N[^-\x)^iX) iV'(-^(A)((.(A)y 


where 


9{T'{u{z)+tV±d)\T') 9{T'{-uiz)+tV±d)\T') 
9{t'{u{z) ± d)\T') ’ 9{t'{—u{z) ±d)\T') 

Notice that from (3.10), for A G CF^\[—M,M] chosen from compact subsets, 




g±27riT'iU«(A) _ g^p 


zfcrt 1 + 0 


Im I In a\ 
In^ a 


u(X) 


, tj 4,0 


and for the same values of A we have at the same time 

1 


u(A) = O ( j — ) = m(oo) 


( 6 . 1 ) 


( 6 . 2 ) 


since with (3.23), 


z{k) = 


\/2 I Int 


Inllncrl 

1 + - -+0 

In cr 


1 

In a 


tr 4 , 0 . 
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Along the line segment [—M, M], by direct computation 

M_(A)=C>^j—O'I 0, X G — r]U [m + r, M), 0<r<- fixed, 

and hence all together 

g±27rir'ty«(A) ^ (T i 0, AeCP^D(0,f), 

for any fixed f : m < f < M, provided we ensure that 

T^=o(l), criO. (6.3) 

In cr 

We shall in fact impose a stronger condition on k than implied by (6.3): suppose subsequently that t = {—x)i 
and i; = — In (l — |sip) are sufficiently large such that 

^a/ 2 - y < >f < ^a/ 2 0<(T<y, /i > 0. (6.4) 

Subject to this constraint we obtain directly with (3.10), 

<±)(A) = l + o(^y^^) =iV'(±)(A), (6.5) 

uniformly for A S CP^\J. Moreover 


Proposition 6.1. As x ^ —oo, |si| f 1 subject to (6.4) we have 

N{X) = g-2,rir'tI/«(oo)<T3Y(^)g2^ir'tI/«(A)^3 = + T(A) (6.6) 

uniformly for A chosen from 911(0, f) for any f : m < f < M. Here, T(A) equals 




where 


i(A) = 


is defined and analytic on C\ [ — 




a(A) + a ^(A) —i(a(A) — a ^(A))A _i2l 

i(a(A)-a-i(A)) a(A)+a-i(A) 

— y 1, A — y oo 


(6.7) 


At this point let us get back to the nonlinear steepest descent analysis. We start from the initial RHP 
1.4 and use the transformation sequence 


y(A) hA A(A) hA Z{X) hA T(A) hA 5'(A) 


just like it occurred in Section 2 and parts of Section 3 in the analysis for >r S [i5, |a/ 2 — <5]. However, instead 
of using the “opening of lens” transformation (3.21) we will not split the jump contours near (—m,m) but 
instead introduce a new model problem on the (contracting) line segment (—r, f) C (—Af, M) with 


In more detail 


r = m 


d > 0. 

fs 


( 6 . 8 ) 


Riemann-Hilbert Problem 6.2. Find a 2 x 2 matrix-valued piecewise analytic function J(A) = J(A;xr) 
such that 

• J(X) is analytic for A € C\(—f,r) 

• The limiting values of J{X) are related via the jump condition 

=-i(^-n(A)) 


J+(A) = J_(A) 


le 


A G (-f,r) 


where k— n(A) G C^{—M,M) has been explicitly computed in Proposition S.j- Also, J±(A) are 
bounded on (—f,f). 
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• As X ^ —oo, |si| t 1 subject to (6.4), 

J(A) = (/ + o(l))T(A) 

uniformly for A G dD(f),f) and T(A) has been introduced in (6.7). 
We choose a solution to this model problem in the following form, 

0 1 

X—fj, 


J(A) = S„(A)e‘T"'^3 / ^ 2^i J-M^ ^ \-u 



A G 13(0,7) : 3A > 0 , , 

-10/ ^ ^ (6.9) 


e i 


A G 13(0, r) : 3A < 0, 


where Bu{X) is locally analytic, 

B„(A) = T(A)e‘i^‘^^ 



e-if<^<^3^ A G Z3(0,r) : 3A > 0 


A G 13(0, r) : 3A < 0. 

Let us quickly verify that (6.9) indeed satisfies the properties of RHP 6.2. First, analyticity of the multiplier 
Bu{X) follows from the jump of T(A), 


T+(A) = T_(A)e‘5' 


0 1 

-1 0 


e-f-3^ AG(-r,r)c (-4,4 


72’V 2 

Thus the jump behavior of J(A) is a consequence of the Plemelj-Sokhotskii formula, 

1 -e-d^-n(A))\ / 0 lA _i 

)(-l o'^ 


7 1 

Secondly, as x —>■ —oo, |si| f 1 with (6.4), 

I'M e-dx-n(M)) ^ 

J-M X- 27ri 27r 


= — In 
27r 


g-t(x-n(A)) jg 
ie 0 


X — m 
X + m 

X — m 
X + m 


rM ^-t(x-n+)) 

m A - /r 27ri 




1 


, A G (-r,r). 

l-M X- 2tt\ 
ln5 t\ 


te In 3 t 


O 


t 2 


uniformly for A G 9Z3(0, r). Here we used a Laplace-type argument for the two remaining integrals based on 
the behavior 


>c-n(A) - 

16 GG- 

- 

1 

(M^ — TO^) 

1 (A — m)7 

A j, TO 

>f-n(A) - 


(M^ — TO^) 

1 (—A — m)5 

, A t —w 


as well as the expansions 


\/2| ln( 


1 -f 


In I Incrl 


21 n( 


O 


1 

IntT 


M= 4 ( 1- 


Thus 


J(A) =[l + 0 


72 

1 


72 | In cr| 

T(A) 


O 


(7 In I In CTI 
In^ a 


+ e In 3 

as t —>■ - 1 - 00 , |si| I 1 subject to (6.4) uniformly for A G dD{0,r) with f as in ( 6 . 8 ). 
Instead of the ratio transformation (3.53), introduce for (7.1), 


CT j, 0. 

( 6 . 10 ) 


i?„(A) = 5(A) 


( 6 . 11 ) 


(J(A))_\ |A|<r 

{P{X))\ \X-M\<r 

(Q(A))-\ |A + M|<r 
^(7V(A))-\ |A| >r,|ATM| >r 

where 0 < r = m + dt~i < | and 0 < r < j remains fixed. The model functions N{X),P{X) and Q{X) are 
(as before) given in (3.26), (3.47) and (3.51). 
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72 



75 


Figure 20. The jump contour of the ratio function i?„(A) defined in (6.11). The two 
vertical lines represent the branch points A = ±m. 


The RHP for i?„(A) is formulated on the contour shown in Figure 20. 
i?(A) the only difference in the jump behavior occurs inside the circle Co 


Compared to the ratio problem for 
centered at the origin, 


{ Ru { X ))^ 


iRu{x))_ 


/ +J-(A) 


/igg3gitn(A) _ ^ 


-iesi(l + - 1 + 


(J_(A)) ' 


for A € (—f, f) as well as 

(R,(A))^= (i?„(A))_J(A) (J(A))-\ AG72niA(0,f) 

(R„(A))^ = (i?„(A))_J(A) J) {JiX)y\ A G 75 n iA(0, f), 

and on the circle boundary dD{0,f), 

(i?„(A))^ = (i?„(A))_ J(A)(iV(A))-\ |A| = f. 

Notice that (since r = 0{t~i) —>■ 0) 

“«=■" (nDii)=(i^)' 

subject to (7.1) as well as 

_ 1 

pm-\-dt 3 

A G (to + dt~3, M — f) : t{K — n(A)) >8t a/(M^ — — m?) d/i > ca/ mt —>■ + 00 . 

J m 

The latter combined with (6.6),(6.8),(6.10), using also the boundedness of iV_(A) away from the three disks 
Dj centered at A = 0 and A = ±M and recalling at the same time our estimations for the jumps given in 
Section 3.9, we obtain 

I|G'r„(-;^: |si|) —t — >■ + 00 , |si| t 1 : ^7/2 — ^ < < -a/2. 


( 6 . 12 ) 
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Opposed to the this L°°-estimation we can obtain better L^-estimations through 



1 _ 

l-f 



2d f"* / 

dA = — (l — cos(27rtT^)) + 2 / (l — cos(tO(A))) dA = O ( 
ts J-m \ 


1 3 In^ t 


and 


dD{0,f) 


|u(A)|^|dA|<^=0 

In t 


1 3 In^ t 


ln3 t 


The last estimation is required for the L^-estimation of the jump on dD{0,f): 

J(A)fV-i(A) = J(A)(T(A))-'T(A)fV-i(A) ^ ("j + o T(A)(fV(A))-' 

V \tB In^ t/J 

^ = ^ (l + 0( ^ \ (l + o( 

V \tB\n^tJJ \ Vt2 In 

= (/ - 27rir'<Fu(A)T(A)a3(T(A))”^ + ...) (^I + O ■ 

We summarize, 

\\GR^i--,t,\si\) - < 1^1 , t-t+oo, |si|tl: “^^2 - ^ < K < ‘^V2, /i > 0 (6.13) 

ts ln 2 t 6 t C) 

and combining (6.12), (6.13) this implies from general theory [12] the unique solvability of the i?„-RHP in 
L^(Sij^), moreover its unique solution satishes 

|i(i?„)_(-;t,|si|)-/|U3(s,j<^^. 

fe ln 2 t 

In the same way as we derived (3.61), we obtain 


r(a:|s) = V—a 


2ti 


^^ 3 ( 0 |r) 6 » 2 (tP|T) ■^^y , ^ ^ 


6 * 2 ( 0 |r) 6 » 3 (tt/|r) 


with 


£ui>() = - [ {Gr^{w) - l)^^dw+ 0 


However, by direct computation. 


/ 1 A i 

r I rl 

Vt 3 Int/ T 

f + / 

^dD{0,f) J — r 


(Gk„ (w) - l)^^dw + O ( — 


1 


O Int 


and thus 


with 


/ {GRjw)-l)^^dw = o(-^) = I {Gr^{w)-I)^^ dw 

J—f \t^ intj JdD{0,f) 

/ 1 X I ^ 1 —k ,,3,^^^/l —kA 1 — kA ^ , . 

■u(a;|s) = -eW-- , , cd 2(-a;)2y K —— ) , —— ) + Ji(a::,s), 


2 vm? 


1 + k/’ 1+k 


(6.14) 


Proposition 6.3. For any given /i > 0 there exist positive constants to = to(fi),vo = uo(/i) and c = c(/i) 
such that 

|-^i(a:, s)| < yt>to,v>vo-. -/i < u < ^^21. 

Together with Propositions 4.2 and 5.1, Proposition 6.3 completes the proof of Theorem 1.10. In order to 
derive (1.29) in Corollary 1.11, we choose t > tQ,v > vq such that 0 < cr < |. Since 


03(O|t) 02itV\r) 03(O|r') O^tVFlF) 


OaitVlr) 6»2(0|r) BoitVr'lF) 6>4(0|r') 

6 » 3 ( 0 |r' 


(6.15) 


04(O|rO 


= 1 + 0 


Ini 


= 1 + 0 


1 


y/i\n^ t 


we get with (3.10) at once 
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Similarly, 


Moving ahead (recall (3.6)), 


e^{tVT'\T') = 1 + 0 


t 


= 6»4(iV^r'|r')- 


^ ^ X 1-k 6i3(0|r) 92{tV\T) _ 

2 Vl + k2 6»3(tM|T) 6>2(0 |t) 




2 \ (—a:)4 ln 2 |a;| 


(6.16) 


uniformly as a; —>■ —oo,|si| f 1 such that 0 < cr < i. In view of Proposition 6.3, the latter estimation 
completes the proof of Corollary 1.11. 


Remark 6.4. The introduction of J{X) in (6.9) together with bounds on N-{X) in ( 6 . 5), ( 6 . 6 ) allowed us to 
control the error Ji{x,s) in the region (6.4). The analysis of Ji(x,s) inside the region 


|(t,u) : t>to,v>vo: y < cr < 5, 


h>0, 0 <6 <- 



which is not covered by Propositions 4-2, 5.1 and 6.3 is much more difficult. In this case (6.3) is violated in 
general and thus N- (A) may become unbounded. For instance, if we were to choose 

1 In < 

t>to, v>vo: -<(J< X— <0, X > 0 , 
then from (3.6) and (3.10), using again (6.15), 


X 

—e-i / —T 


1-k 6»3(0 |t) 02{tV\T) 


2 v/m? 03 {tV\T) 02 (O|t) 


2 




Hence we only restore the leading Hastings-McLeod asymptotics for x < | • For the values o/ x > | more 
terms will contribute to the leading order, or equivalently, the leading asymptotics changes once the “Stokes 
line ” 

2 /A 1 , 

V = - v2 1 — - Int 

3 6 

is crossed. We provide further detail on this interesting feature in Section 8 below. 


7. Regular transition analysis near and above the separating line xr = 1-^/2 
In this section we assume that both t = {—x)i and v = — ln(l — |sip) are sufficiently large such that 

k>‘^V 2-^, f2GR. (7.1) 

Again, we use the transformation sequence 

y(A) hA X{X) hA T(A) 

but subsequently the nonlinear steepest descent analysis of the T-RHP in (2.2) is very different from the 
one used in Sections 3 and 5. In fact, the steps carried out below are a modification of the ones used in the 
analysis for |si| = 1, see [19], chapter 11, §4. The details are as follows. 

7.1. g-function transformation. We hx 

in (2.2) and define for A G C\X1t with St 
tj7(A) = T{X)e^<^sW-^W)<y3.^ ^ith 



= [-A*, A*] U 72 U 75 U 72 U 75 U 7fc, 

5(A) = I (a 2 - 0 ' = d(A) + O (A-i) , A^oo. (7.2) 
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This transforms the T-RHP to a RHP with jumps on [—A*, 0] U [0, A*] such that 


'Ct7_|_(A) = t37_(A) 
ri7+(A) = ■n7_(A) 


(1 - Si + si(l - S 1 S 3 ) 

-S 3 (1 - SiS3)e“™(A)y ’ 

(1 - siS 3 )e™(^) -S3 

Si + si(l - S 1 S 3 ) (1 - siS 3 )e 


A e [-A*,0] 
-tn(A)^ I Ae[0, A*] 


where 


n(A) = ff+(A)- 5 _(A) = 


3 V 2 


-A' 


>0; A e - 


1 1 


For the (ll)-entries in (7.3) and (7.4) we use (3.1) and observe that 


>!r-n(A)>0 |A| > Ao = 




V2 


\ 2 V 2 J 


(7.3) 

(7.4) 


Hence the strict inequality for xr — n(A) will fail on the (shrinking) interval [—Ao,Ao]. We thus face the 
necessity to consider an additional model problem at the origin. Besides this new feature, we however still 
expect (and will prove) that the outer model problem is given by a problem in which we neglect the diagonal 
entries in (7.3),(7.4) and apply estimations (3.2) and (3.3) for the off-diagonal, i.e. on the full line segment 
[—A*, A*] we require a function T(A) such that 

TRA) = T_(A,(» “). A.(-TF). (7.5) 

Furthermore, along the remaining infinite branches, i.e. along 

= 72 U 75 U 72 U 75 U y 7 fc, 
the jumps in the S'-RHP are once more of the form 

ti7+(A) = tI7_(A)e-‘9(^)'"3GT(A)e*9(^)'^^ A S 

Hence by triangularity and the sign of 3?(g(A)) these jumps are asymptotically close to the unit matrix away 
from the endpoints A = ±A*. We proceed now with the analysis of the relevant model problems. 


7.2. The outer parametrix. The outer problem consists in finding a 2 x 2 piecewise analytic function 
T(A) = T(A; e) such that 

• T(A) is analytic for A G C\[— 

• If the cut is oriented from — ^ to the boundary values of T(A) are related via 

rRA).T_(A)(» “), Ac(-^,y (7,6) 

• T(A) is square integrable on the closed interval ;^] 

• As A —)■ 00 , 

T(A) = /-tC>(A-i) 

A solution to this quasi-permutation RHP appeared already in Proposition 6.1, formula (6.7). Note also that 

--I-O (A“^) , A 00 . 


T(A) =/ 


(7.7) 
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7.3. Origin parametrix. Our construction is motived by the local expansion 

>f-n(A) = >r-|v^+2v^A2-^+C>(A®), A ^ 0, 

3 2v 2 

and solves the following problem. 

Riemann-Hilbert Problem 7.1. Determine a 2 x 2 pieeewise analytie function F{X) such that 
• F{X) is analytic inside the disk D{0,r)\{—r,r) for some fixed 0 < 


(7.8) 


Orienting the segment (—r, r) from left to right, we have the jump condition 
F+{X) = F_{X) 


g-t(>t-n(A)) jg 


, AG(-r,r); n(A) = ^ Q - A^ 


ie 0 

• As X ^ —oo, |si| t 1 subject to (7.1) we have 

FiX) = iI + oil))T{X) 

uniformly for |A| = r > 0. Here, T(A) is the outer model function introduced in Section 7.2. 

The difference between this problem and RHP 6.2 is only marginal: Recall that cr = — k and set 

0 l' 

F^^{C) = ^ 2i4 /r 7^'^^ 


(7.9) 



e <i 


-I 4 <^<^3 

SaC < 0 . 


(7.10) 


We require the following Lemma. 
Lemma 7.2. Ts C —>■ oo, we have 


m = 


27ri 


o-u 


d/4 = 


u-C 1 -xHL 




G 


Proof. Assume first that C ^ hence 


m = - 


27riC , 




^clC = - 


1 

27riC 


/ d/i + 0 (C-i) 

.JR 




Secondly for ^ G K we apply Plemelj formula 


4(0=±i«-‘’ + ^v.p,/^ 


c 


d/i. 


The principal value integral can be simplified as follows: let denote the contour shown in Figure 21 with 
0 < e < i?. 


A — .R + i 


X-R 




A + R “t" i 


A — £ A + 6 


A + R 


Figure 21. Integration contour y/j ^ for a possible choice A > 0. 
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From Cauchy’s theorem 


On the other hand 


Vi? > e > 0 : 


7 -R, 


c 


d/4 = 0 . 


o-^J■ 


c 


d/i = —i / e 




as well as 


and 


K+fl+i rC-R \ g-M" 




-R+i / M ~ C 


d /4 = i [\-iC+R+it?^ 


hm 

R—^oo 

eiO 


/O 


r-c-e /-C+^A 


Thus 


4(C) = T 




C--R ic+e J C 
_^2 J 1, 3C > 0 _ i/tt 


i? + Jq 

dfi = v.p. 


+ i r e-(C-^+i^)^ 

Jo 


dt 


R—>oo 


R — it 




M- c 


d/4. 


:d /4 + e i ’ ’ ' “ = ——— (l + O (C 0) + s 

^ ' 0, 3C < 0 27riC ^ 


5C > 0 

1 0, 5C < 0 


27ri Vjj+i /4 - C ^ 

which completes the proof. 

The solution to the model problem is now given by 

F(A) = So(A)F«^(C(A)), |A| <r; C(A) = V? (a + >r - n(A))^ = (2^2 (^1 - ^ - ^A" + O (A®)^ 

where we introduce the (locally analytic) multiplier 

Bo{X) = T(A)e’J-3 I 0 ) " ^ ^ : 5A > 0 . ^ 1 -k 

A e iA(0,r) : SA < 0 ^ ^ 

More detail on the local behavior of Bo{X) near A = 0 will be needed later on, first 


□ 


Bo{X) = 


j_(i ‘ 


72 VI 1 y 2 v-1 1 

31i 
M 


Vl 1 


A" + 


8 V-1 1 


Ov + |^(_\ ;)v + 0(A«) A^O, 


64 Vl 


followed by 


{Bowy 


2 , 3i z' 1 — 1 


A" + - 


1 1 


Vi(-i O+sG i 

T#(-i 


and thus 

5o(A)cr+(i?o(A)) ^ 2(1 _jg^ “ ^'^ 2 A +-eusA^ - ^cr 2 A^ + —ecrsA'^-^cr 2 A® + C> (A®) . 

We directly verify the required jump condition (7.9), and as t —>■ 00 (subject to (7.1)), 

Fix) = (i- iee"‘/(C(A))iio(A)(T+(iio(A))-') T(A) 

= (/+^^i?o(A)(T+V(A)+0(t-i))T(A)= (i + o(t-5))T(A), a+= (7.11) 

uniformly for 0 <ri<|A|<r 2 <;^ with rj > 0 fixed. 
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7.4. Edge point parametrices near A = ±A*. The construction is almost identical to the one presented 
in Section 3.8. More precisely for the right edge point we take 

A"(A) = |A-A*|<r, (7.12) 

where {() was introduced in (3.45) and we have 

B, 3 (A) = T(A)e*^-3('^. -Ij «(A)-3 ^ |A-A*|<r 

with the locally conformal change of variables 

C(A) = (^|e-t 5 (A)) ' ~ V2i2t)l (a - , |A - A*| < r. 

The reader easily checks that A’'(A) defined in (7.12) has the same jump behavior as shown in Figure 12 
with M replaced by A* and g{X) as in (7.2). Moreover 


A’-(A) = 


/ + T(A) 


48C 


i /-I 6e 
-6e 


^/)+0(C-S)|(T(A))-^ 


T(A), C^oo 


so that for t —>■ +oo, |si| f 1 such that k > 

A"(A) = (/ + 0(t-i))T(A) 

uniformly for 0 < ri < |A — A*| < r 2 < 

Near the remaining edge point, we choose 

A^(A) =B4(A)1^^(C(A)) + 

7=1 


with Al(C) as in (3.49) and 




= T(A)e'J-3 


, |A + A*|<r 


where we use the change of variable 


3t 


(7.13) 


(7.14) 


C(A) — (^-^g(A)j ~ A/2(2t)3 ^A + ^j , |A + A*|<r. 

The jump behavior of (7.14) is depicted in Figure 14 with —M replaced by —A* and the g-function as in 
(7.2). We note that 


A^(A) = 


/ + T(A) 


1/1 6e 


; +o(c-m(T(A))-^ 


48/5 V“6e - 
hence, as t —>■ +oo, |si| f 1 such that > \\/2 — 

A^(A) = {l + 0{t-^))T{X) 


T(A), C^oo, 


(7.15) 


uniformly for 0 < ri < |A + A*| < r 2 < This completes the parametrix construction and we can move 
ahead with the ratio transformation. 
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7.5. Final transformation - ratio problem. Define 


C(A) = w{X) < 


(A’'(A))-\ 

l^-^l 

< r 

(A^(A))-\ 

l^+ 

< f 

{F{X))-\ 

|A| < r 


(T(A)) \ 

-H 

> r 


(7.16) 


where 0 < r < ^ and 0 < r < | remain fixed throughout. This leads to the RHP shown in the Figure 
below. 


72 



Figure 22. The jump contour of the ratio function ^(A) defined in (7.16). 


• ^(A) is analytic for A G C\E5 where Ej = [—A*, A*] U Co U Cr^ U Ci^ U Es^ 

• The jumps are described as follows: Inside the two circles resp. we simply copy the cor¬ 
responding expressions from Section 3.9, namely the jumps inside Cr 2 resp. Ci^ with the proper 
identification of the gr-function (7.2) and the replacement n{X) = 0. The same applies to the infinite 
branches 72 U 75 outside Cq as well as for 71 U 73 U 74 U 75, with the replacement N{X) i-)- T(A) as 
determined through (7.6). On the line segments (—A* -I- f, —r) U (r, A* — f), observe that 


^+(A) = e-(A)T_(A) 

^+(A) = e-(A)T_(A) 


_ige-*(^-nO))\ 

ies3 J 

_igg-t(><-n(A)) \ 


(T_(A)) \ AG (-A*+f,-r), 
(T_(A))-\ AG(r,A*-f) 


and inside Co, first for A G (—r, 0), secondly for A G (0, r), 

e+(A) - e-(A)F_(A) ( _.^^_t(.+n(A)) i, ,, ^ ) (^-(A)) 


-1 


e+(A)=e-(A)F_(A) 


le S3 


ies3 -f 

(—ie — 53)6“*^^“^^''')^ 
-lee _|_ e 

Notice that for A G (—A* -I- f, —r) U (r, A* — r), 

1 


(J^-(A)) 


-1 


t(>r - n(A)) > t{>c - n(±r)) > “/2 + t ( ^7/2 - g ( 2 “ 


> -/2 + ds,t, ds> 0 
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and for A S (—r, r), 


< dge 


-|V2t 


(si(l + e-^*) - k)e-‘("-n(^)) < dioe-s^* 

But since T_(A) is bounded on the full segment [—A*, A*] and i^-(A) inside (—r,r), we obtain that 

(7.17) 


|si|) ~-f||L 2 nL“([-A*,A'']) ^'^iie t —>■ c», |si| t 1 : /2 € 


The jumps along 72 U 75 inside Co are again exponentially close to the unit matrix, hence we now 
list the contributions on the circle boundaries and Cq, 

^+(A) = e-(A)A’'(A)(T(A))-\ A e Cr,-, i+{X) = ^_(A)A'(A)(T(A))-\ A e 

and 

$+(A) = e-(A)F(A)(T(A))-\ AS Co. 

From (7.13),(7.15) with fixed radius f. 


AeC,3UC^3: C^A) = / + C (t-i) , <^oo,|si|tl: he 


On the other hand with (7.11) and again fixed radius r, 


A € Co : Cj(A) — I -\- O 2 ^^ t—>'CX),|si|j'l: 3 ^^ -^2 € K- 


(7.18) 

(7.19) 


• As A —)■ c» we have that ^(A) —>■ I. 
We solve the RHP for ^(A) iteratively. 


7.6. Iterative solution. The following estimation is obtained as in Section 3.9, we only have to keep in 
mind (7.17),(7.18) and (7.19), 

IIQ(-;i, |si|) - ^l|L 2 nL“(S£) < ^+ 00 , |si| t 1 : > ^a /2 - y, /2 e K- 

Thus the ratio problem has a unique asymptotic solution in L^(Sj) which satisfies 

IIC-(sii Isil) - < dt~^, d > 0 . 

7.7. Extraction of asymptotics - proof of Theorem 1.12. Through the transformation sequence 

y(A) hA A(A) hA Z{\) hA T{\) hA w{X) hA ^(A) 

we obtain 

u{x\s) = 2 lim fA(y(A;x)),J = 2^=^ lim \zhz)T(z)e-*^3{z)-H^))a3 

A->oo z-yoc L J 12 

Since 

g(z)-d(z) = ^+0(z-^), T(z)=/-^ui+C(z- 2 ), z 

we deduce from the integral equation for ^(A) the exact identity 


00 


i(a;|s) = V—3 




c(r 


where (as t —?> 00 , |si| t 1 such that >f > | a /2 — ^,f 2 e K), 

^(^) = {G^{w)-l)^^Aw + 0{t-^) = 

and the last contour integral was computed by residue theorem, compare (7.11). Hence, as a: —>■ — 00 , |si| f !> 


r(x|s) = -e\ -- 1 - 


(7.20) 


with 
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Proposition 7.3. For any given /2 € K there exist positive to = to{f 2 ),vo = vo{f 2 ) and c = c(/ 2 ) such that 


\j2{x, s)| < ct ^ yt>to, V >Vo, V > -V^t — /2. 

o 

This completes the proof of Theorem 1.12 and at the same time allows us to derive Corollary 1.2. 


7.8. Proof of Corollary 1.2. Suppose both t and v are sufficiently large such that 

.. In(l-kiP) . 2 

- } Is > o V 2 . 

( —x) 2 o 

In this case a < 0 is bounded away from zero and hence (7.19) can be improved to 

G^(A) =1 + 0 (f-se'"*) =1 + 0 (r°°), A e Co, 

compare (7.11). Therefore, 


|si|) ~ — Tj t ^ +00, |si| t ^ ^ 


and (7.7) gets replaced by 


£{>c) = 


t 



t3 


+ o {r ^). 


with X and 7 independent coefficients Uk- These could in principle be computed by residue theorem in (7.7) 
using (7.13) and (7.15), but it is easier to work with the differential equation (1.10) itself: As a: —>■ — 00 , |si| f 
1 , 


u(ai|s) 



(7.21) 


uniformly for x > Ia/S and the asymptotic expansion can be differentiated with respect to 


X. 


Remark 7.4. The modulus |si| will appear explicitly in (7.21), but only contained in exponentially small 
error terms. All power like terms are \si\-independent and in fact identical to the ones derived in [13], 
Theorem 1.28. 


We go back to the Fredholm determinant and note that for fixed 7 < 1, 


^ Indet (/ - jKAi) 


L^{x,oo) 

Hence with Remark 1.6 and (7.21) after integration 


pOO 

/ Mis(y;7)d2/= «s(a:;7))^ - a;uis(a;;7) - Mis(a;;7), x S K. 

J X 


Indet (/ — 'yKAi) 




(7.22) 


where 0 ( 7 ) is x independent and there exist universal to,vo > 0 and d > 0 such that 

k(x, 7 )| < Vt = (-x)5 > to, -ln(l- 7 ) > ^ 0 , x>^V 2 . 


Now keep t >to fixed and let 7 f 1 in (7.22), so that with (1.7), 

limc( 7 ) = Inco = In 2 + C'(-l). 
7ti 24 


This, together with (7.22), completes the proof of Corollary 1.2. 
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8. Appearance of “Stokes-lines” - proof of Theorem 1.13 
We assume that t = {—x)i and v = — ln(l — |sip) are sufficiently large such that 


x> /jSK. 


Working with this scale, several steps in the nonlinear steepest descent analysis of Section 7 can still be 
applied. In fact, the only difference appears in the matching relation (7.11). Subject to (8.1), we have that 

< max {l, } 

and thus the two model functions F{X) and T(A) are in general no longer asymptotically close to each other 
on dD{0,r). In order to deal with this feature, we choose G Ig = (—oo, |((7 + 1)) C M with q G Z>o and 
notice that for C ^ K, 


9-1 


/(c) = ^/^dM = -^Ec 

27ri 7 e M - C 27riC ^ 


n—O 

L^J 


f d/r + ^ 

TO 27ri , 




n—0 


’■-q 


2ti\ 


C 


C 


d/x 


d^. 


which implies with Lemma 7.2, as t —>■ 00 subject to (8.1), 

[ g-i j 

e'^‘/(C(A))=-gJ ^ C-^-i(A)r(i+n) +0(t 

n=0 ^ ' 


max{/3-i-L2±iJ).-i-L2^J} 


and the last term is small because of our assumption /a G Iq. We now factorize (7.11), 
F(A) = (/ - iee"‘/(C(A))Bo(A)i7+(Bo(A))-') T(A) 

= Bo(A) E r"”-'(A)rQ+n^a+j (i?o(A))-' 


( 8 . 2 ) 


xi?o(A) U-^e'^V«(A) 


dMa+I (i?o(A)) 't(A) ^ (A(A) x F(A))T(A) 


m-C(A) 

and go back to Section 7. After the explicit transformations 

r(A) hA A(A) hA T(A) hA w(A) hA ^(A) 
we use another: With (7.16) and Figure 22, let 

oa3_/e(AW), |A|<r 
\^(A), |A|>r 

so that by definition, the function 0(A) solves the following (singular) RHP. 

Riemann-Hilbert Problem 8.1. Determine a 2 x 2 matrix valued function 0(A) such that 

• 0(A) is analytic for A G <C\{Yj^ U {0}) 

• Along we have jumps, 0a(A) = 0_(A)G'e(A), with 


\g^{X), XGCe,UCr,UX:, 


(8.3) 


• 0(A) has a pole of order 2[2^J +1 at X = 0 provided q G Z>i. In order to be more precise about 
the singular structure, notice that for q G Z>i, 


e(A) = 0(A)(F;(A))-' = 


"( 9 ) ^ "( 9 )-! 

E (' 

e=o 


k=l 


in(q)+e) 


X^ 


A=o ("(l) + ^)! 


(EW) 


-1 
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with n{q) = 2[2^J + 1, is analytic at X = 0. Thus with (8.2) and comparison, we can express the 
full singular part {0_fc} in terms of the first n{q) terms of the regular part of 0(A). This is done 
explicitly for q € {1, 2 }, 

^(7t 


0_i = res 0(A) = -(A0(A))' 


1 



a=o2 


-ie) 


V 


p = 


(8.4) 


s/2 271 V 2 (ii/2)2 

but becomes much more involved in case q G Z> 3 ; we shall postpone the general case q G Z >3 to a 
forthcoming publication. 

• As X ^ oo, we have 0(A) -G I. 

The jump matrices in RHP 8.1 are asymptotically close to the unit matrix, however 0(A) has a pole at 
the origin. Such type of problems have been analyzed in nonlinear steepest descent literature, see e.g. [8, 9], 
at least for the case of a first order pole. Following loc. cit., we first note that (8.4) is equivalently stated as 


0(A) = 0(A) 


1 -a]t- 
0 1 ' 


|A| < r; 


T= — (^ 

72 V-ie 1 


K = 


1 _ IP. 
V 2 


where 0(A) is analytic at A = 0. Hence, 0(A) is indeed characterized uniquely by the properties stated 
in RHP 8.1. Secondly we can resolve the singular structure for q G Z>i via another transformation, the 
dressing transformation, 

0(A) = (a^L^J+7 + + ... + aHi + Ho) S(A) ^ , A e C\S 5 (8.5) 




21 I 

where are A-independent. The function S(A) is characterized as follows. 

Riemann-Hilbert Problem 8.2. Determine a 2 x 2 matrix valued function S(A) subject to the following 
conditions 

• S(A) is analytic for X G C\I]j 

• We have that S+(A) = S_(A)G'e(A) for X G Ej, with Gq{X) given in (8.3). 

• Ts A —>■ oo, 

S(A) = / + 0(A-i). 

Notice that the unknowns Bj are algebraically determined from the singular structure in the 0-RHP and 
(8.5), i.e. for g S {1, 2} we have 

-1 


Bo = -^m 


le 1 
1 -ie 


:(0) + -S'(0) 


ie 1 
1 -ie 


or for more general q G Z >3 equations for the remaining unknowns in term of ^(0) and higher derivatives 
thereof. Observe that all jumps in the S-RHP are close to the identity since 

\\Gei-;t, |si|) - /|U2nL~(E,) < c > 0 

uniformly for t —)• oo, |si| f 1 such that x: > 172 — /s^ and —oo < f 3 < ^{q + 1). 

8.1. Extraction of asymptotics. The RHP 8.2 is iteratively solvable and we have for its solution 

kil) -< ct— 

Notice that from a direct computation, by tracing back the transformations, 

F(A) hA X{X) hA Z{X) hA T(A) hA tJ7(A) ha ^A) ha 0(A) ha S(A), 

we get 


u(a7|s) = \/—x 
atinr 

w(a;|s) = y/—x 


2H7„-i 


+ -/ (^-{w){Ge{w) - 1 )) dw 


qGZ 


>0 


■^2 

We will now continue with the special case q G {1,2}, i.e. we compute 

^ 1 ^ ^^max{ 2 / 3 - 3 ,/ 3 -f ,-3}^ 


H-2 = 0. 


72 


2 i?7- 


1672t 872^2 


( 8 . 6 ) 
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Notice that 

and thus 


:(0) =1 + 0 (^tmax{/3-|,-i}^ ^ 2'(0) = O -1}^ 

1 'j ^^max{/3-§,-l}^ 


Bn = -- 


(8.7) 


From this it follows in (8.6) that 


M(a;|s) = —£\—-z 


2 \1-^^ 


ep 


1 Itt e° 


+ J3ix,s), - = - — 


^ — 1 / 1 , cr = - \/2 — K 

27rV2(2v/2t)2 3 


and 


Proposition 8.3. For any given /a € (—oo, g) there exist positive to = io(/ 3 )j'^’o = ^oifs) and c = c(/ 3 ) 
such that 

2 2 

|>^ 3 (a;) s)| < ct“ yt>to, V > vq, V > -V^t — fslnt. 

o 

If we were to choose /a G (—oo, g), then p = 0(tB-i^ = o(l) and we restore the leading order Hastings- 
McLeod asymptotics, 

m(x|s) = I Ja(a:, s)| < 


In other words, 


2 r- 1 

Si: V = —v2t -Int 

3 6 


forms the first Stokes line and we have completed the proof of Theorem 1.13. 

8.2. Proof of Corollary 1.17. Observe that (8.6) and (8.7) allow us to derive a slightly better expansion 
than the one given in Theorem 1.13, as a; —>■ —oo, |si| f !> 


u(a;|s) = -e \/-2 




where, uniformly for xr > |v^ - /a^, /s G (-oo, ^), 

I I Jmax{t-2e2'^‘,t-3}, /g g (-oo, i] 

Next, we use identity (1.4) and combine it with (8.8), 


( 8 . 8 ) 


(8.9) 


02 

-^Indet (/-yi^Ai) 


L‘^{x,<x>] 


= WAs(a;;7) = -ohi—5F + 


8.^ 


X —>■ —oo 


and from (8.9), noting that e°’‘ = e with v fixed, 

Jj |i7i(a;,7)|d2a; = O ^ 

However, by straightforward computation. 


, kiltl, >f>^v^-/3^, /sG^-oo, 0. (8.10) 


dx'^ 




and i? 2 (x, 7 ) satisfies also (8.10). Hence, we have shown that for x —>■ —oo, |si| f Ij 

_ x^ 1 

L2(a:,oo) 12 8 


Indet (/ — yTfAi) 


ln|x|+aa: + & + ln(^l- + C> ^ 
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uniformly for >c > /s G (—00, |) with some ^-independent constants a = 0(7), b = These 

can be computed as in the proof of Corollary 1.2, i.e. we use 

d 2 

— lndet{l-jKAi) = 7)dy = «s(a^; 7)) - xul^ix;"/) - xGR. 

OX L^(ai,oo) 

together with (8.8) to conclude 0(7) = 0 and the limit 7 t Ij with t > to fixed, yields 

lim6(7) =lnco = In 2 C'(-l)- 

7ti 24 

9. Singular transition analysis for k g [<5,1-^/2 — 5] with 0 < <5 < fixed 

Several steps of the analysis carried out below have their natural counterpart in the regular analysis. One 
major difference however lies in the construction of a new outer parametrix which naturally involves the 
exceptional set as defined in (9.7) below. First we let 

>f=^S(0,oo); t = {—x)^ > 0, 0 < !> = — In (|sip — 1 ) = —27r/3, 1 < |si| < 

and keep >c G [(5, — <5] with 0 < (5 < hxed. Also here, the following expansions will prove useful. 


9.1. Preliminary expansions. Note that 




|si| = l+-e--‘+0(e 


— 2>ct \ 


and 

Si = e”’ 

Moreover, 


(^1 + ‘ + O = S3, Si + si(l - S1S3) = (^1 - + O . 

arg(si) = f (1 + f e—‘ + O (e-^-*)) 


so that 

Proposition 9.1. As t —>■ 00 , |si| j, 1 with e‘ 2 *^ = ie, 


Si — ie ^1 - (1 -|- ies2) e + O (e , S3 — —ie 2 ~ 16^2) e + O (e 


Si -h si(l - S1S3) = ie ( 1 - - (1 - ies2) e 


0 \ 


uniformly for x G [6,\^/2 — 5] with 0 < (5 < fixed. 


9.2. g-function transformation. The bijective correspondence between x G (0, fv^) and k G (0,1) has 
already been established in Proposition 3.2, and all local expansions in Corollary 3.3 can directly be applied 
to the singular transition analysis. With this in mind we define (formally just as in (3.12)), 


g{z)=4i - m^))^ dfj., zGC\[-M,M] 

Jm 


where 


1 


1 


1 


k 1 , .r 1 

0 < m = —= . ^ < - < M = —= , ^ 

V 2 Vl + k2 2 Vl + k2 ^2 


and k = ^{x) is determined through (1.22). The relevant analytical properties of g{z) are summarized in 
Proposition 3.4 and we choose in the T-RHP (2.2) (as before). 


After that, introduce, 


A* = M. 

S{X) = ^ A G C\Et 
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where we shall use the same notation S{X) for the transformed function as in the regular transition analysis. 
Observe that, for A G {m,M), 


g-t(^-n(A)) 


S+{X) = S.iX) 
so that (as before) 

Gs{X) 

Similarly, for A G (—M, —m) 
S+iX) = S.{X) 
hence (as before) 

Gs{X) 


-S3 

_ p-xt\ _p-t(xr+n(A)) 


si(l —e ^‘) —e 


'j , >f — n(A) = 8 / ■\/ (M2 — — m?) d/j, > 0, 

/ J m 


0 

0 


—)■/, a; —)■ — cxD, |si| 1 : kG 


5,-V2-6 


S>0. 


—e n(A)) 


_ggg-itn(A) 




_g— t(>t+n(A)) 


)\ . - 

J , K — n(A) = 8 y \/— /a2)(^2 _ 


0 

_g—i^e—itn(A) 


gif e+itn(A) 

0 


/, X —>■ —c», |si| ), 1 : >f G 


5,-y2-(5 


, (5 > 0 . 


The only difference occurs in the gap, first for A G (0, in) 

-1 


_g3gitn(A) 


5+(A) = S.{X) 

- - siS3))e“‘*^(^) l) ( 


(si + Si(l - SiS3))e -(siS3-1)2 

-1 0\ fl S3e“(^)\ 


0 -1 


0 




= 5_(A)e‘^"^5L,(A)e-‘^"^5B(A)e*^'^^5c/,(A)e-*^'^^ AG(0,m) 

and secondly for A G (—to,0), 

c QN _ c -1 (si+si(l - SiS3))e“(^)' 

f.+ (A) - i,_(A) -(S1S3-1)" 


= ()( 


0\ f-1 0 

0 -1 


1 -(si + Si(l - siS 3 ))e“^(^)\ 

0 1 


= 5'_(A)e‘^'^3S'L,(A)e-‘^'^^5'^5(A)e*^'"3S'c/,(A)e-*^‘"^ AG(-m,0). 

Here, 5 'lj(A) and Suj{X) have appeared in (3.18), (3.19) and Sd{X) = —I. Hence, we can use Proposition 
3.6 in the given context as well and open lens: Compare Figure 5 and define 


L(A) = 


This leads to the following RHP 

Riemann-Hilbert Problem 9.2. Determine a 2 x 2 matrix valued function L{X) such that 

• L{X) is analytic for X G C\([—M, —m] U [—m, m] U [m, M] U U 7 ^ U 

• The jumps are as follows, see Figure 5 for orientations, 


5(A)e*i'^3 5-i(A)e-if-3^ 

XGCt 

5(A)e*5"35^^(A)e-if-3^ 

XGCf 

5(A)e'i'^3^-i(A)e-if-3, 

X G £2 

5(A)e'i"3^^^(A)e-if-3^ 

A G £2 

5(A), 

otherwise. 


Gs(A), 

X G (—M, —m 

e*t‘"3 5'^^(A)e-if'^3^ 

Ag 7 + 

e*t‘"3 5'^^(A)e-if'^3^ 

A G 73 " 

e*5‘"3 5'j^^(A)e-if'^3^ 

AG72+ 

e*t‘"35'^^(A)e-if‘^3^ 

AG72” 

Sd{X), 

A G (—m, m) 


L+{X) = L_{X){ 


(9.1) 
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• A —>■ OO, 

L{X) = I + O {X-^) . 

The jumps in this problem along the lens boundaries and along the infinite contours are exponentially 
close to the unit matrix in the (singular) double scaling limit as long as we stay away from the branch points 
A = ±TO, ±M and A = 0. Compared to the L-RHP in the regular transition analysis, the outer parametrix 
will now also have a jump in the gap (—m, m). 


9.3. New outer parametrix. In case of the singular transition analysis the outer model problem is posed 
on the full segment {—M, M) and reads as follows 

Riemann-Hilbert Problem 9.3. Determine a 2 x 2 matrix-valued piecewise analytic function N(X) = 
N{X; t, e) such that 

• N{X) is analytic for X G C\[—M,M] 

• We have jumps 

0 

Yje—itn(A) Q 

A G (—m, to) 

_g-i5 eN 

0 


N+{X) 

= p- 

V-e 

N+{X) 

= 7V_(A)e*’"‘"^ 

N+{X) 

= »-(A)(/5. 


gif e+itf2(A) 


, A G (—M, —to) 


, A G (to, M) 


(9.2) 

(9.3) 

(9.4) 


• N{X) has boundary values on [M, M] 

• As X ^ OO, 


N{X) = I + O {X-^) 

Notice that the function 

N{X) = A G C\[-M, M] 

has jumps on {—M, M) described as follows 

pitn(A) 

0 


iV+(A) = iV_(A) ( 


A G (—M, —to) 


iV+(A) = iV_(A)e*’"‘"^ A G (-TO, to) 

0 1 


iV+(A) = 7V_(A) 


-1 0 


A G (to, M). 


Up to the presence of the diagonal jump in the gap (—to, to), the jumps are thus identical to the ones which 
we generated in Proposition 3.11. To compensate for the diagonal part, we slightly modify our construction 
of Section 3.5: Let u{z) denote again the Abel (type) map 

tt : CP^\[—M, M] —>■ C, z i-G-u{z) = j tv = u{oo) - \-O (^z~^') , z —>■ oo; ^ ~ ~ 


IM 


for which we summarized certain properties in Proposition 3.9. We require 

. (y-wNAw), (..) 

where (formally as before, see also (1.23)), 

V{k) = V = —— [ (M^ — — p?) dp = —n(z), zG{—M,—m); d 

J-m 27r 


r 

'4 


and 


d(z\T) = 03 (z|t) = ^ exp [ink^T + 27rifcz] , z G C 


fcGZ 


is the third Jacobi theta function. With this, we have the following analogue to Proposition 3.11 for the 
current situation (9.2)-(9.4). 
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Proposition 9.4. The function 

]\T(\) = „-i7rM(oo)a-3 f W4>W ^2 gi7ru(A)(T3 -ifetra ('9 6') 

e{^+tv) 

is single-valued and analytic in CF^\[—M,M]. Its jumps are stated in (9.2),(9.3) and (9.4), furthermore, as 
A —>■ oo, 


iV(A)=/+- 


e'a+tv) 

— C-KT^-ri^ — ITTC 


_ ^(0) 0(u{oo)— ^ —tV—d) —27riu(oo)gi^e 


_ e( 5 +iy) ^ e{^+tv) e(«(oo)-d) 2 i 

\ I ^(0) 0{u{oo)-\-2-\-tV—d) M—m 27riufooI„ —iHe ( 2 “^^^) 1 {-rn 

' e(j+ty) e(ti(oo)-(;) 21 ® ® '-e(f+ty) “t i^rc 


+ 0(A-"), 

provided (x, si) are uniformly bounded away from the discrete set 


Zr,. = 


|(a;,Si): 2tV(x) = n G Z\{0}'j. 


(9.7) 


Proof The stated jump behavior follows from a direct computation using Proposition 3.9 and the properties 
of Jacobi theta functions. For the normalization, we take into account the shift by the half period appearing 
in the numerators in (9.5), hence we have to guarantee that 

which is ensured for the values {x, si) away from the set Z„ given in (9.7) as the simple roots of 6{z) = 93 {z\t) 
are located at z = A + | mod Z + tZ. □ 


Remark 9 . 5 . The latter Proposition is in sharp contrast to Proposition 3.11. There we do not require the 
half-period shift and thus we always have that 9(tV) 0 since tV{>c) G K, i.e. there is no exceptional set. 

From now on we shall always assume that (x, si) is bounded away from the set Zn. 

9.4. Local parametrices. The required model functions near A = 0 and A = ±m slightly differ from the 
ones used in the case of the regular transition. Near A = 0 we analyze 

Riemann-Hilbert Problem 9 . 6 . Determine a 2 x 2 piecewise analytic function H (A) such that 

• 11(A) is analytic for A G (C n D(0,r)^\(^j^ U U (—r,r)) withO<r<Y- 

• For the boundary values, 

H+{X) = i?-(A)(^J j, AGD(0,r)n(j+Uj+) 

H+(A) = J), AeiA(0,r)n(7rU72-) 

H+(A) = -H_(A), Ae(-r,r) 

• As t ^ 00 , |si| 4 , 1 such that k G [J, — <5], uniformly for A G dD(0,r), 

II(A)=(l + o(l))N(A). 

A solution to the latter problem is derived in two steps. First define 

argCG (0,f)U(2f,7r) 
argCG (-f,0)U(7r,f) 

^ j, argCG(f,2f) 

1 ). «sC€ (-¥,-!) 


^°(C) = 
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with r2(0) = 7 rV(>t:) from (1.23) and where we orient all six rays from zero to infinity. Second, the required 
parametrix itself: With JV(A) from (9.6), 


H(A) = N(A) 


+1, S5A > 0 
-1, S5A < 0 


X i^°(C(A))e-‘^(^)‘"^ AeD(0,r), 


(9.8) 


where 


t 


((A) = - {n(A) - fl(0)) = StMmA {l + O (A^)) , A G L>(0, r). 


The required jump behavior follows from a direct computation using a local contour deformation, moreover 
we have 

HiA) = {l + 0{t-°°))N{A) 

as t —>■ +00, |si| j. 1 such that k G [<5, fv^ — 5] uniformly for 0 < ri < |A| < r 2 < y. 


Near A = m we have to modify (3.33). Consider 

'l, 


'1 1 ' 

,0 1 




;:). 

iW 1 0^ 

0 ij i-1 1, 


argCG (0,f) 
argC G (%,7r) 


argC G (tt, 


argCG (f,27r) 


with Alo(C) as in (3.34). This leads to the following bare RHP 

• is analytic for C G C\{argC = 0, ^,7r, 

• We have jumps on the contours sketched in Figure 23, 


H«^(C) = H«^(C)(^J argC: 


2?! 

Y 


• As C —> 00 , 
A«^(C) = < 
Now define 

where (as before) 


H«"(C) = AY(C) 
Hj^(C) = 


A^^iC) = A^^iC) 


+1, 3C>0 i 

-1, 3C < 0 ^ 2 Vl -i 


-1 0 
0 -1 


, arg C = TT 


1 0 

1 1 

-1 1 

-1 0 


dir 

, argC=y 


, argC = 0 


48CJ 


-1 -6i 

-6i 1 


O 


(o 


17(A) = H,3(A)A^^(C(A))e-5«"(^)'^^e-'^''^^ A G Z?(m,r) 


efC^-3 


(9.9) 


C(A) = 


6te‘5 / ((M^-^2^ (m^-^2^)5 d^ 


2 

^ (A — to), |A —m 


— m\ < r 


Br,{A)=N{A) 



3A > 0 
S5A < 0 


X e 






A-M \ 
A — m J 


6 {A) = 


/ A — TO A 

\a-m) 


and 
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argC = f 



Figure 23. Jump 
contour for the bare 
parametrix 


0 1 



,-t(»-n{A)) 


Figure 24. Jumps of the model 
function U{X) near A = m. 


is analytic near X = m. This leads to the jump behavior of U{X) as indicated in Figure 24, moreover we 
have the matching relation, 


C7(A) = 


I + N{X) 


1 f-1 6e 


48CS 


+0(rS) (iV(A))-^ 


NiX) 


and thus, as t —>■ +oo, |si| j, 1 such that >c S [J, — J] with 0 < J fixed, 

UiX) = {1 + 0(1)) NiX) 


uniformly for 0 < ri < |A — to| < r 2 < min{^, ^{M — m)}. 
For the remaining endpoint A = — m, we put 

'-L 


24""(C) 


1 0 

^g-i7r(l-7) Ij ’ 

^ 1 o\ 

^g-i7r(l-7) j ’ 

^ 1 o\ A 

g-i7r(l-7) ijlo 


1 


argC e (-7r,-f) 
argC e (-f,0) 

argC e (0, f) 

argC e (f ,7r) 


with Alo(C) as in (3.41). We obtain that 

• A^^iO is analytic for ( S C\{argC = —1,0, |,7r} 


• Along the curves shown in Figure 25, 

A^^iC) = A^^{0 {„-mi- 

\RH 


1 


e *..( 1 - 7 ) 1 


, argC = -- 


RH ,/ / — 1 0 


A«"(C) = A""(C)( argC = 0 




RH { 


Af^iO = _,-i.a-7) 


1 _gi’r(l-7) 

0 1 
-1 


, argC = 


gi7r(l-7) 

0 


, argC = 7r 
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ai-gC= i 




Figure 25. Jump 
contour for the bare 
parametrix 


Figure 26. Jumps of the model 
function V (A) near A = —m. 


As C —> 00 , 


1-1, 9C<0 2 Vl 


+ o(c 


efiC 2 <T 3 ^ 


jgi7r(l-7) 


48C5 


1 


For the parametrix itself, we put 

F(A) = B£3(A)I^"(C(A))e-3'‘^®(^)‘"^e-'T^<^3^ g D[-m,r) 


(9.10) 


with 


C(A) = 


6t {{M^ - d^i 


2 

^4t\/2m(M2 — ^ (A + to), |A + m| 


< r 


and the locally analytic multiplier 

'+1, 3A>0 


B,,iX)=NiX) 


-1, S5A < 0 


X e "1 


F(A) = 


/ + N{X) 


— ie‘’^(l-T) _igi’r(l-7) 

-1 1 


0ggi7r(l-7) 


J(A)-'^3 


,X + M\~^ 
X + m 


48Ci 

so that as t —t 00 , |si| ), 1 such that xr S [J, — <5], 

F(A) = (J + o(l))7V(A) 


';-^'Uo(ri) (iv(A))-‘ 


»(A) 


uniformly for 0 < ri < \X + m\ < r 2 < minl^, |(M — m)}. This completes the construction of the relevant 
local model functions. For A close to ±M, we can simply work with the same functions P{X) and Q{X) as 
in (3.47) and (3.51). 


9.5. Ratio transformation. We work with the ratio function i?(A) defined in (3.53) but with H{X) from 
(9.8), U{X) from (9.9), F(A) from (9.10) and A^(A) as in (9.6). By construction, R{X) is characterized by the 
following properties. 

• R{X) is analytic for A S C\E/j where the contour is shown in Figure 15. This follows directly 
from our construction of N{X) and H{X) as both functions precisely match the jump behavior of 
L{X) in the gap (—m,m), compare (9.1). 
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• The jumps of R{X) are only slightly different from the ones listed in Section 3.9. We choose not to 
list all expressions explicitly, instead only collect the important estimation 


R+iX) = R.iX)GR{X-,t,\si\), 


l|GR(-;t, |si|) - IWL^nL'^i^R) < p c > 0 


(9.11) 


as t —)■ + 00,1 Si I I 1 such that xr € [<5, — 5] for 0 < J fixed, provided we stay away from 

the singular points (cc, si) G Zn- 
• As A —>■ oo, we have that R{X) -g I. 

With the help of estimation (9.11) the RHP for i?(A) can be solved asymptotically in L?'{T,r) and we have 
for its unique solution 

| si |) - i\\l^(^r) < 

9.6. Singular transition asymptotics. Just as in the regular case, we have the exact identity 


i{x\s)=2V^ lim 


= Z\/—xe 


^2ti 




where 


^ J (^R_{w){Gr{w) - l)^ = O [t ^) . 


Notice that with (3.24) and identities of theta functions, 

(M i = 6>(0) 0{u{oo) -^-tv -d)M- TO 27 rm(oo) iHe ^ ^ 3 ( 0 |t) 6 > 3 (tH|r) M - TO ■ 

^ ^ 0{^+tV) 0{u{oo) — d) 2 i 6 * 2 ( 0 |t) 02 (tH|T) 2 i 

and hence we obtain 


i.(a:|s) = 


„2U 


—e(M — m) 


6 ' 3 ( 0 | t ) e^{tV\T) ■^y 


02(O|t) 6>2(tH|r) 

Now using Proposition 3.16, combining it with (1.34) and the identities 




(9.12) 


M — m = 


1 1 -k 

vi+k2 ’ 


6*2(0|t) 

6»3(0|t) 


1 -k 
1 + k’ 


one deduces 


and 


M(a;|s) = — ^ ^ ^ dc ( 2tVK 


2 \/m? 


l-kA 1 -k 


1 +k/’ 1 + k 


+ J4(a;,s), 


Proposition 9.7. For any given 6 G (0, Iv^) there exist positive constants to = to{S),vo = vo{S) and 
c = c(S) such that 

\Ji{x, s)| < ct“s, 'd t>to, V > Vo '■ tS < V < t ^ 2^/2 — 
provided (a;, si) are uniformly hounded away from the exceptional set Zn introduced in (9.7). 


10. Extension at the lower end of singular transition 

The relaxation of the lower bound for x is handled via the same approach as in Section 4, i.e. for t = (—x) 5 
and V = — ln(|si|2 — 1) sufficiently large such that 

2 /- 

0 <t ’*<>c<-v 2 — J, with fixed 0 < ly < 1 , 
o 

we obtain with the help of contracting radii the following analogue of (4.4) 

' 1 -kA 1 -k 

Vl + kj ’ 1 +k 


m(x|s) = —C\j — -R , dc (2{—x)^VK 


+ 0({-x) 


\-i+U 


( 10 . 1 ) 
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uniformly as x —^ —oo, |si| 1 such that 0 < for any 0 < ly < | away from the exceptional 

set (9.7). After that we would use the analogue of (4.6), respectively Remark 4.1, for the singular analysis 
and this combined with improved estimations gives 

£(>f) = C> +c>(t-i). 

We summarize. 

Proposition 10.1. For any given 5 G (0, |\/2),i7 G (0,1) there exist positive constants to = to{6,r]),vo = 
vo(S, rf) and c = c(i5, if) such that 

v2-6y 

provided (x, Si) are bounded away from the exceptional set (9.7). 

In order to obtain estimation (1.38) in Corollary 1.15, we would choose t > to such that 0 < xr < t~^ and 
again (x, Si) stay away from (9.7). Then, by Corollary 3.3 we have 


J4(x, s) < Ct 


'ifl.f-i';} vt>to, v>vo- t^~'^<v<t\^ 


X 1 + k / 3 /l-k\ 1 -k 

eW-- , ^ dc 2 {-x) 2 VK 


2 Vl+W 


1+k/’1+k 


cos W: 


tV{>c)) 


O 


((-X) 10 ) 


Next, for 0 < XT < t = with (3.10), (compare (4.7)), 

TTtV{>c) = —^(—x)5 — /31n ^8(—x)5^ + /31n I In (|sip — l) | — /3(1 + ln27r) + O 
and from Stirling’s formula 


argT 



argT 



/31n I In (|sip - l) I - /3(1 + ln27r) + O (t 


( 10 . 2 ) 


as 


xt = 0{t^) —>■ + 00 . Together 


TrtV{>c) = —-(—x )2 — /31n ^8(—x)^^ — ip 




(-) 


where (compare Theorem 1.20), 


+ = -arg T ( 2 + ) - arg si. 


We substitute (10.3) into (10.2), use the addition theorem for cosine as well as the identity 

+ ^^ = —esin^;, zGC 


cos 


to obtain, as x —>■ —oo, 


(10.3) 


—e< / — 


dc(2(-x)5+A: 


2 vTTl? 


l-k\ 1-k 


1 + k/ ’ 1 + k 


/—X 


sin ^|(—x )2 + /31n (8(—x) 2 ) + +0 {{—x) lo^ 

+ o((-x)-ra), (10.4) 


4 

uniformly for D < x < t~'s and away from the zeros of the trigonometric function appearing in the last 
denominator. Hence, the Jacobi elliptic function leading term reproduces also here the known singular 
oscillatory leading order of (1.20) for (t,x) such that t > to and 0 < x: < t~s. Next we will also obtain 
control over the error J 4 (x, s) in the same region. 
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11. Further extension at the lower end for singular transition 

We go back to Section 5.3 and assume throughout that both t and v = — ln(|sip — 1) are sufficiently large 
such that 

t>u'=+^>0, t>to, fcGZ>3. (11.1) 

The difference between regular and singular transition asymptotics can be read off from (5.3), 




( 11 . 2 ) 


Hence we can still use most of the nonlinear steepest analysis worked out in Section 5, however a crucial 
differences occurs in the matching relations (5.4) and (5.6), compare [8] where the same approach was first 
used in the derivation of (1.20) for fixed |si| > 1. We shall factorize (5.4), 


4''-(A) ^ Er{X) 


/ + Cr_(^(A)) X! (('^+ l)2m - -hcr+(^(A)) 




m—1 

^(A)-2m+l 


m! 2" 


E 


{v)2m 0 

0 (_)™(_j,_l)2^ 


C(A) 


— 2m 


ml 2" 


4-^ (A) 


with 


/ 1 0 

ErW = 

V *3 Cl-^) 


1 


This factorization is necessary since with (11.2) 

= ('c(A)-^ + ^‘ 


C(A) 


A-A'* 


— 2vq 


A +A* 
A-A'* 




i.e. the multiplier Er{X) is not close to the unit matrix. Similarly, 

ml 2 

A —2m+l 


^^(A) ~ Ei{X) 


u _\A —2m—1 

I + Cr+ (/3(-A)) ^ E - (u)2m)--^ (/^(“A)) 

m—1 


2^ 




m—1 


m—1 




C(-A) 


— 2m 


ml 2” 


^^(A) 


with 


Ei{X) = a2Er{-X)a2. 

This means that the ratio problem (5.7) for x(A) cannot be solved iteratively at this point and we have to 
use further transformations. First, compare [8], equation (38), we let 

rx(A)ii;.(A), |A-A*|<r 
‘I’“(A) = E(A)F;KA), |A + A*|<r 
[x(A), |ATA*|>r 

where we use (this is in fact the only difference to [8]) a shrinking radius 

1 1 L 

- > r = - > t fc+i. 

4 V 

The Riemann-Hilbert problem for $“(A) was analyzed in [8], its jump contour is identically as shown 
in Figure 19 but with additional first order pole singularities at A = ±A*. We choose not to reproduce the 
analysis here, compare equations (40) — (42) in [8]. Subsequently the singular structure was resolved by 
another transformation 


$“(A) = (A/ + H)$‘^(A) 


(A-A*)-i 0 

0 (A + A*)-!'’ 


A G c\j:^ 
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where the A independent matrix B is determined algebraically as in equation (46) in [8]. It is well defined 
for all values (cc, Si) bounded away from the discrete set 


-t-\- j3 ln(8t) -\- (f = nTT, 

O 


(11.3) 


compare again [8], shortly before Remark 3. If we agree to stay away from the latter exceptional points, the 
jumps for <I>^(A) have the following structure 


4>i(A) = <l>l(A)f(^ ^ 


»'! —1 


0 


-1 / A-A* 

0 


0 

A +A‘ 


, A e 7j, j = 


and 


$t(A) = <1>1(A)((^ 


-1 X - X* 


-1 x - x * 


$t(A) = <i>l(A)(^(^ ^^P^,^_^Ei\X)¥{X){^^{X)) 

Hence we have 

G^.(A)=/ + o(e"^-‘-'=‘^'), AG7,, J = 
so that (just as before in (5.8)), 


A + A^ 


A +A^ 


XeCr 

XeCi. 


For X e Cr we estimate 

|G^.(A;t,|si|)-/| <c 

with 


G^d(-;t, |si|) — /||L 2 nT“(U 7 j) ^ d^e 

{xtf 


fc-1 
— dai ^ + 1 


k G Z>3 


(11.4) 


rt 






e* e 


xt 


{xty 


S{X-t, |si|) 


0{r) = — arg C(A) 


A + A* 

'a-A* 


Thus, using again bounds for the error term £ from [29], we obtain that with the contracting radius 

||G$d(-;t, jsij) —/||i2pij;,cx,(c,,uCi) ^ *=+i, t > toi fc G Z>3. (11-5) 

Combining estimations (11.4) and (11.5) we can thus solve the ih^^-RHP asymptotically subject to the scale 
(11.1) as long as we stay away from the discrete set of points (a;, si) defined via (11.3). The derivation of the 
leading order follows now as in [8], we only have to adjust the error terms according to (11.5): from equation 

(51) in [8], 


i(a:|s) = 2 y—xBi 2 + O ^(—a:) ^(fe+i) ^ ^ k G Z> 4 , 


where we used (11.5) and the a-priori estimation for the unique solution of the <I>‘^-RHP. The matrix entry 
Bi 2 is then computed explicitly in [8] and we only have to adjust the error terms. Hence, 


r(a:|s) = 


sin ^|(—a :)2 -|- /3 In ^8(—a:) 2 ^ + + O ^(—a:) 


/ 2fc-7 \ 

+ 0 1{-X) 2()= + l)j 


as X —>■ —oo subject to the scale (11.1) for k G Z >4 provided we stay away from the points (x, si) defined im¬ 
plicitly via (11.3). The latter estimation, with fc = 4, appears in Corollary 1.15, expansion (1.38). Moreover, 
with (10.4), we have now shown that 


m(x|s) = -ej-| -^=^dc ( 2{-xyVK 


l-k\ 1-k 


1 -fk / ’ 1 -Fk 


-I- J4{x,s), 


and 


Proposition 11.1. There exist positive constants tg and c such that 

I J4(x, s)| < yt>to, 0<v<t^ 

provided (x, si) are uniformly hounded away from the set (9.7). 
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This concludes the proof of estimation (1.36) in Theorem 1.14, we just have to combine Propositions 10.1 
(with f? = |) and 11.1. 


12. Extension at the upper end for singular transition 
We choose t>to and v = — ln(|sip — 1) sufficiently large such that 

^^2- ^ < K <\y/2 ^ 0<a=\V2- /i > 0. (12.1) 

o c o o t 

This implies with (3.10), 

0 < |2tE(>!r)| < c > 0 universal. 

In other words, we can always guarantee that for sufficiently large t and — ln(|sip — 1) subject to (12.1) the 
exceptional set Zn is empty, and thus u(a;|s) is asymptotically pole-free. This is important as we can now 
go back to Section 6 and simply repeat our line of reasoning, modulo minor adjustments. First, from (9.6), 

AT(\\ = ^(Ok') -2TTiT'u(oa)tVa^ ( 2mT'u{X)tV<73 


valid for any A € C\[—M, M] and where 




0(t'(u(z) +tV+^± d)\T') 0(t'(-m(z) +tV+^± d)\T') 


9{t'{u{z) ± d)\T') 


6(t'{—u{z) ± d)|r') 


But 


e(T' (u{z)+tv +'^±d^ |r') =e^{T' {u{z) + tv ±d)\T'), z e CF^\[-M,M], 
and thus with (6.1) and (6.2), subject to (12.1), 


= 1 + 0 


Incrl 


= N',^^\X). 


Summarizing, Proposition 6.1 also applies to the singular transition x —>■ —oo, |si| j, 1 subject to (12.1). Now 
following the logic of Section 6, we would employ the sequence of transformations 


y(A) hA X(A) hA Z(A) hA T(A) hA 5'(A) 

but not open lens, compare (9.1). Instead, just as in Section 6, a new parametrix J(A) is introduced on 
the contracting segment (—f, f) C {—M,M) with f as in (6.8). The difference between J(A) and J(A) as 
considered in RHP 6.2 is minor, we would simply enforce 

J+(A) = J_(A)(^ y, Ae(-f,f). 

This difference amounts to an additional sign in the (12)-entry of (6.9) which contains the Cauchy transform. 
However, subsequently all steps are identical to Section 6 and we obtain 

“(xl.) = (2(-x)*rif (1^) , 1^) + J4(x,.) 

with 


Proposition 12.1. For any given /i > 0 there exist positive constants to = to(/i)j vq = 'Co(/i) c = c(/i) 
such that 

|J4(x,5)| < Vt > to, /4 <V< ^V2t. 

' ' In t 3 3 

This completes the proof of Theorem 1.14 and expansion (1.39) in Corollary 1.15 is derived as in the 
regular transition case. 
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13. Singular transition analysis near and above the separating line = Iv^ 

In order to derive Theorem 1.16 we would choose identical steps as those in Section 7: We work with 
the same outer parametrix T(A) as in (7.6) and the same endpoint parametrices A’'(A) in (7.12) and A^(A) 
in (7.14). The only (minor) difference occurs in the definition of F{X), i.e. in the definition of the origin 
parametrix. For the singular transition we require the jump behavior 


F+(A) 


F_(A) 


_„-i(x-n(A)) 


le 



A e (-r,r); 


n(A) = 


8 

3 



3 

2 


instead of (7.9). This difference is resolved by an additional sign in the (12)-entry containing the Cauchy 
transform in (7.10). After that we have to replace the matching (7.11) through 

FW =[l- + O (f^)) T(A) 

and then simply copy all subsequent steps carried out in Section 7. This leads us to 

/ I X I ^ f 1 T / 2 /- / x3 

u(a;|s) =-e^-- [^1 + —+ J 5 (a:,s)j , a=-V2-K, t={-x)^ 


with 

Proposition 13.1. For any given /2 € K there exist positive tg = ^o(/ 2 ))L'o = Vo{f 2 ) and c = c(/ 2 ) such 
that 

\j 5 {x,s)\ < Ct~^ Vt>to, U > Ug, V>‘^V2t-f2, 
and thus completes the proof of Theorem 1.16. 


Appendix A. Expansions of complete elliptig integrals 
We frequently use the elliptic integrals £’(k) and K(k) defined as 

dt 




0 V(l-k2t2)(l_i2) 


, ke (0,1), k' = x/i - kA 


i-F 

For k ), 0 or k t 1, the expansions of £'(k), K{k) are well known, compare for instance [29]: 

^ ln(l - k^) + dim)) , 0<x/r^<l; 

m\m\ V 2 / 

m=0 ^ 2 


1 “ (1) iSj 


with 


E{k) = 1 + - V ^ ^ 

2 (2)mw! 

m=0 


(1 - k^)’"+^ — ln(l - k") + d{m) - 


1 


(2to + l)(2m + 2 )) ’ 


d{m) = -0(1 + to) - 0 - + m , d(0) = 2 In 2; (a)n = 


F(a + n) 

r(a) ’ 


a 7 ^ 0,-1,-2,... 


and 0(ai) denoting the digamma function. Furthermore, 

OO /lx /■ 1 




K(k) = - y 

2 to ! to ! 

m—0 


r7.n.x ’’’ ( 2)"i(2)™i_2m 

^(k) = 2 to! to! 

m—0 


|k|<i. 
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Appendix B. Jacobi theta and elliptic functions used in the text 


The Jacobi theta functions 6 * 2 (z); ^ 3(^)1 ^ 4 (^) used in the text are defined by the formulae (see e.g. [29]) 
where > 0 , 

00 

6^{z\r) = 0{z\t) = 1 + 2 ^ cos(27rA:z), 

k^l 

00 

^(_l)fcgirfrcos(27rfcz), 

k^l 

00 

= 2 ^ gi7r(fc+i)V + l) 7 rz) . 

k=0 

The zeros of the theta functions are located at 



r =1 + 2 


6»4(z|r) = 93 [ z +- 


92{z\t) = + I 


up to shifts by vectors of the lattice A = Z + rZ; all zeros are simple. Now set q = and denote 
9j{z\T) = 9j{z,q). With 


q = exp 


K'jk) 

A(k) J 


k = 


92(0, q) 
93(0, q) 


e (0,1), 


the Jacobi elliptic functions cd(z, k) and dc(z, k) are given by the formulae 

^' 3 ( 0 , 9 ) 02(C, 9 ) 1 . 2 


cd(z,k) = 


c = 


02 (O,g) 03 (C, 9 ) dc(z,k)’ " 2 A(k) 

for all values z S C away from the zeros of @ 3 , resp. 02 appearing in the denominators. Finally we note that 
by an application of a Landen transformation, 


, K{k) _ K'{\) 
^A'(k) K{\) ' 


A=^, A' = +-A2. 
1 + k 
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